﻿ANALIZA CIRCUITELOR ELECTRICE CU PSPICE LUCRĂRI DE LABORATOR/ CONSTANTINESCU FLORIN, IONESCU ANISOARA, MARIN CONSTANTIN-VIOREL, NITESCU MIRUNA CUPRINS Prefața 1 Breviar teoretic 3 2 Utilizarea programului Pspice 33 3 Lucrări de laborator 61 4 Soluțiile problemelor 89 5 Bibliografie 132 CAP 1 BREVIAR TEORETIC I CIRCUITE REZISTIVE 1 Existența și unicitatea soluțiilor 1 1 Circuite liniare Se spune că un circuit rezistiv are soluție unică dacă ecuațiile acestuia sunt satisfăcute simultan de o mulțime unică de tensiuni și curenți Există circuite simple care nu au soluție sau au E1 - E2 un număr infinit de soluții Circuitul din figura a are o soluție unică I = ——— dacă RxO Dacă R=0 se pot distinge doua situații: I a b i) dacă ErE2=0 rezultă 0*1=0 și circuitul are o infinitate de soluții (orice Ie R ) ii) dacă ErE2xO rezultă 0*1x0 și circuitul nu are soluție Circuitul din figura b are o soluție unică U = *s1 + *s2 dacă GxO Dacă G=0 se pot distinge două G situații: i) dacă Isi+Is2=0 rezultă 0*U=0 și circuitul are o infinitate de soluții (orice Ue R ) ii) dacă Isl+Is2x0 rearltă 0*Ux0 și circuitul nu are soluție Condiția de existentă și unicitate a soluției unui sistem liniar de n ecuații algebrice Ax=B este detAd) Această condiție exprimă faptul că ecuațiile sistemului sunt liniar independente intre ele Dacă detA=0 există doua situații de interes practic: i) rangul lui A este n-l și determinantul minorului principal bordat cu coloana termenilor liberi este nul; în acest caz una dintre ecuații este liniar dependentă de celelalte și sistemul are o infinitate de soluții, ii) rangul lui A este n-l și determinantul minorului principal bordat cu coloana termenilor liberi este nenul; în acest caz sistemul nu are soluție în circuitele simple prezentate mai inainte se observă: i) în figura a, pentru R=0 avem detA=0 și există o bucla formată numai din surse de tensiune 3 ii) în figura b, pentru G=0 avem detA=O și există o secțiune formată numai din surse de curent Aceste observații pot fi generalizate pentru orice circuit format din rezistoare dipolare cu R>0 și surse independente, condițiile de existență și unicitate fiind exprimate în funcție de parametrii circuitului Teorema Un circuit liniar format din rezistoare dipolare liniare cu R>0 și surse independente are o soluție unică pentru orice valori ale tensiunilor electromotoare și ale curenților electromotori dacă și numai dacă sunt satisfăcute următoarele condiții: - nu există nici o buclă formată numai din surse de tensiune, - nu există nici o secțiune formată numai din surse de curent Demonstrație: Necesitatea se demonstrează prin reducere la absurd Intr-o buclă de surse de tensiune din teorema a II-a a lui Kirchhoff rezultă £Ek = 0 Dacă Ek sunt astfel încât £Ek 0 circuitul nu are soluție Dacă Ek sunt astfel încât EE] 0 Rezistorul strict crescător are Ли^Лі>0 Rezistorul de tip U (de la “unboundcd” - nemărginit) arc o caracteristică cu proprietățile: -dacă este controlat în curent atunci lim u(i) = ±~ i^±~ -dacă este controlat în tensiune atunci lim i(u) = ±~ u^±^ Rezistorul de tip H (de la “half-unboundcd” - scmi-ncmărginit) arc o caracteristică cu proprietatea rezistorului de tip U numai pentru una dintre extremități (spre sau spre -«>) Vom enunța în continuare, fără a le demonstra, o teoremă de existență ,o teoremă de existență și unicitate și o teoremă de unicitate 6 Teorema 1 (de existența) Fie un circuit format din surse independente și rezistoare dipolare pasive controlate în tensiune sau în curent astfel încât variabila de control poate lua orice valoare între și +~ Acest circuit are un număr impar de soluții dacă sunt satisfăcute următoarele condiții: i) orice buclă formată din rezistoare controlate în curent care nu sunt controlate și în tensiune conține cel puțin un rezistor de tip U sau două rezistoare de tip H orientate diferit în raport cu sensul de parcurgere al buclei, ii) orice secțiune formată din rezistoare controlate în tensiune care nu sunt controlate și în curent conține cel puțin un rezistor de tip U sau două rezistoare de tip H orientate diferit în raport cu sensul de parcurgere al secțiunii Observații : i) având un număr impar de soluții circuitul are cel puțin o soluție; de exemplu dioda tunel din circuitul din figura 3 este un rezistor pasiv controlat în tensiune iar circuitul are trei soluții ii) restricția cu privire la bucla formată numai din surse de tensiune din teorema pentru circuitele liniare este un caz particular al condiției i (sursele de tensiune pot fi considerate rezistoare neliniare controlate în curent care nu sunt controlate și în tensiune) iii) restricția cu privire la secțiunea formată numai din surse de curent din teorema pentru circuitele liniare este un caz particular al condiției ii (sursele de curent pot fi considerate rezistoare neliniare controlate în tensiune care nu sunt controlate și în curent) Condițiile de existență și unicitate a soluției unui circuit neliniar cu rezistoare crescătoare pot fi formulate cu referire la existența și unicitatea unor circuite liniare Teorema 2 (de existență și unicitate) Fie un circuit N format din surse independente și rezistoare crescătoare controlate în tensiune sau în curent astfel încât variabila de control poate lua orice valoare între - și + «> Considerăm circuitul N’ obținut prin înlocuirea rezistoarelor neliniare cu rezistoare liniare cu rezistențele Rk > 0 Dacă N’ are o soluție și numai una pentru orice valori Rk și în N’ există o pereche arbore-coarbore astfel încât rezistoarelor controlate în curent din N le corespund rezistoare din arbore în N’ și rezistoarelor controlate în tensiune din N le corespund rezistoare din coarbore în N’, atunci N are o soluție și numai una Observații i) condițiile de existență și unicitate a soluției privitoare la sursele independente sunt aceleași ca la un circuit liniar; ii) spre deosebire de teorema precedentă, condițiile sunt strict topologice (nu există condiții legate de forma caracteristicilor rezistoarelor neliniare); ca urmare observăm că un circuit cu rezistoare crescătoare se comportă similar cu un circuit cu rezistoare liniare din punct de vedere al 7 existenței și unicității soluției Teorema 3 (de unicitate) Fie un circuit format din rezistoare dipolare crescătoare controlate în tensiune sau în curent astfel încât variabila de control poate lua orice valoare intre și +~ Dacă există, soluția acestui circuit este unică dacă sunt satisfăcute următoarele condiții: i) orice buclă formată din rezistoare controlate în curent conține cel puțin un rezistor controlat în tensiune ii) orice secțiune formată din rezistoare controlate în tensiune conține cel puțin un rezistor controlat în curent Observații : i) restricțiile cu privire la bucla formată numai din surse de tensiune și la secțiunea formată numai din surse de curent sunt cazuri particulare ale condițiilor i) și ii), ii) teorema nu asigură existența soluției, de exemplu dacă în circuitul din fig l a se înlocuiește dioda Zener ideală cu o diodă Zener reală (care este un rezistor controlat în tensiune) atunci acest circuit satisface condițiile teoremei dar nu are soluție, iii) un rezistor strict crescător poate fi considerat controlat atât în tensiune cât și în curent; un circuit format din surse independente și rezistoare strict crescătoare satisface condițiile teoremelor 1, 2 și 3 dacă sunt satisfăcute restricțiile cu privire la bucla formată numai din surse de tensiune și secțiunea formată numai din surse de curent, deci are o soluție și numai una 2 Scrierea ecuațiilor potențialelor nodurilor 2 1 Circuite liniare Fie un circuit care conține rezistoare liniare dipolare, rezistoare liniare multipolare având o reprezentare controlată în tensiune și surse independente de curent Pentru început se consideră circuitul pasivizat Fie i vectorul curenților și fie и vectorul tensiunilor laturilor grafului acestui circuit Toate elementele au o reprezentare controlată în tensiune, deci se poate scrie i, = Gt -u unde Gt este matricea conductanțelor laturilor Dacă A este matricea de incidență laturi-noduri teorema I a lui Kirchhoff se scrie A- i = 0 sau AGfU=0 în circuitul inițial la care se conectează și sursele independente de curent, teorema I a lui Kirchhoff se scrie Ai = ig unde is este un vector ale cărui componente sunt sume ale curenților surselor independente conectate la fiecare nod Rezultă A-G|-u=( și deoarece u=Al-V (V fiind vectorul potențialelor primelor n-l noduri iar Vn=0) rezultă A-GfA’-V =is sau notând matricea conductanțelor nodurilor cu Gn,1=AG1A‘ rezultă Gn 1-V=iscare este ecuația potențialelor nodurilor Atât pentru scrierea ecuațiilor de către un operator uman, cât și pentru construirea matricei 8 Gjj-j într-un program dc calcul, sc pot stabili reguli mult mai simple decât utilizarea formei matriceale prezentate anterior Să considerăm un exemplu simplu pentru care vom serie teorema I a lui Kirchhoff Alegem V3=0 în nodul de potențial Vp (V1 - V2 )G2 + (V1 - V2 )G3 + (V1 - V3)G1 = G(y\ - V2) în nodul de potențial V2: (V2 - V1 )G2 + (V2 - V1 )G3 = +1 Sau VAG, + G2 + G3 -G)- V2(G3 + G2 -G) = 0 W +2G3) +v2(G2 + G3) = +1 V3 = 0 Trecând în membrul drept termenul produs dc sursa comandată ecuațiile devin: V1(G1 + g2 + G3)- v2(g3 + G2) = G(V1 - V2) • - v1(g2 + G3) + V2(G2 + G3) = +1 V3 = 0 Sc poate imediat deduce o regulă simplă dc scriere prin inspecție (dc către un operator uman) a ecuațiilor potențialelor nodurilor Ecuația în nodul j este dc forma: V L Gk Jkk 2 ViGk = k & ,j k &j !sk unde prima sumă conține toate conductanțclc conectate în nodul j, a doua sumă conține toate produsele ViGk unde Gk este conectată între nodurile j și i iar suma din membrul drept conține toți curcnții surselor (independente sau comandate) conectate la nodul j (luați cu semnul + dacă intră în nodul j și cu semnul - dacă ies din acesta) Sc observă că o conductanța G conectată între nodurile i și j apare în matricea G^ în patru poziții cu semnele: coloana/ i ,■ /linie 1 J i + - j - + similar sc poate observa că o conductanță dc transfer a unei surse comandate ca în figură apare în Gn j tot în patru poziții cu semnele din tabel 9 coloana/ / linie i j k + l + Pentru un program de calcul este mult mai simplu să introducă fiecare conductanță proprie sau de transfer în cele patru poziții precizate de aceste reguli decât să construiască matricole A și Gt și să opereze cu acestea Până acum am considerat că circuitul nu conține surse ideale de tensiune (sau, dacă le-a avut, ele au fost transformate în surse ideale de curent) Există însă și surse de tensiune care nu se pot transforma în surse de curent Dacă circuitul conține doar o singură sursă de acest tip atunci se alege Vj =0 și rezultă V2 =E (nu se mai serie teorema I a lui Kirchhoff în nodul de potențial V2) Evident, dacă există mai multe surse de acest tip care au un nod comun sau care formează un “mini-arborc” nodul de potențial nul va aparține acestei structuri iar potențialele celorlalte noduri ale structurii respective se pot exprima ca sume de tensiuni electromotoare Dacă circuitul conține două astfel de surse care nu au un nod comun se alege V^O, rezultă V,=E și se introduce o necunoscută suplimentară / Când se scriu ecuațiile în nodurile de potențiale V3 și V4 latura cu E, se consideră parcursă de curentul / Acestei necunoscute suplimentare îi corespunde ecuația suplimentară V4-V3=E2 Până acum am considerat că avem doar surse comandate în tensiune Comanda în curent se poate transforma în comandă în tensiune: Dacă această transformare nu este posibilă (ca în figura de mai jos) se consideră o sursă cu E= 0 în latura de comandă și se procedează ca în cazul sursei ideale de tensiune (în ecuațiile nodurilor de potențiale V[ și V2 latura de comandă apare prin -i și +z și se adaugă ecuația V^V,) 10 Deci algoritmul de scriere a ecuațiilor potențialelor nodurilor este: • se fac toate transformările posibile ale surselor de tensiune în surse de curent și ale comenzilor în curent în comenzi în tensiune, • se alege potențialul de referință astfel încât cât mai multe potențiale ale nodurilor să poată fi exprimate ca sume dc tensiuni electromotoare, • considerând și necunoscutele suplimentare (curenții unor surse de tensiune conectate între alte noduri decât cclc dc la punctul precedent și curcnți dc comandă) sc scrie sistemul dc ecuații: V] Z Gk k ej k k ei ,j к ej Jsk S VGk = la care se adaugă ecuațiile suplimentare de forma V| - Vm = Ea Exemplu: V3 V4 = 0 V2 V1 = E1 V1 R1 1 1 -1 R2 " V2 V^ = 1 R2 V3 = 2 V1 V2 R1 R4) ' = h 12 După ce s-au determinat potențialele nodurilor sc pot calcula tensiunile ca diferențe ale potențialelor și apoi curenții 2 2 Circuite neliniare Fie un circuit cu rezistoare controlate în tensiune și surse independente dc curent Scriind ecuațiile constitutive ale rezistoarelor i=g(u), teorema I a lui Kirchhoff A 4= 0 și teorema a II-a a lui Kirchhoff u=A‘V rezultă ecuațiile potențialelor nodurilor Ag(A‘V)=is Sursele de tensiune care nu sc pot transforma în surse dc curent și comenzile în curent sc tratează similar cu circuitele liniare 11 3 Determinarea unei soluții prin metoda Newton-Raphson Fie un sistem de ecuații algebrice ncliniarc f(x)=0 unde x este vectorul necunoscutelor Relația care definește metoda iterativă Newton-Raphson se obține dezvoltând pe f în scrie Taylor în jurul lui x® (x la iterația j): unde matricea J(x®)= f ^x(J + 1) j = f A Hf1 dfn dx1 dx1 df1 dfn dxn dxn x=x(j) x ( J + 1) — x + sc numește Jacobianul sistemului și se calculează în punctul x® Impunând condiția ca x!j+i> să fie soluție a ecuației f(x)= 0 și neglijând termenii de ordin superior din dezvoltarea în scrie Taylor sc obține relația de recurență a metodei iterative Newton-Raphson: x(j+1) = x(j) — J—1(x(j))f(x(j)) Algoritmul pleacă de la o aproximație inițială x'"' și efectuează un număr fix de iterații No Dacă eroarea între două iterații succesive x(N) x(N—1) este mai mică decât o eroare impusă £ se consideră că metoda este convergentă iar soluția este x-N (||X| = i n 2 X x ) Dacă după No iterații k=1 k eroarea nu scade sub valoarea £ se consideră că metoda este divergentă Dacă f(x)=0 arc mai multe soluții, prin această metodă sc determină numai una dintre ele și anume cea care este mai “apropiată” de aproximația inițială xfU) Dacă x arc o singură componentă metoda Newton-Raphson arc o interpretare geometrică simplă Fie f(x) având reprezentarea grafică din figura de mai jos și fie xfU) aproximația inițială Valoarea x(1) se obține prin x(1) = x(0) — f1—1 (x(0)) • f (x(0)) adică, deoarece tga = f’(x(0)) = f()x (1) xfl) sc poate determina intersectând tangenta la f(x) în xfU) cu axa Ox 12 în continuare sc obțin xf2), x1'” și soluția numerică tinde către soluția exactă x* Metoda Ncwton-Raphson nu este convergentă pentru orice aproximație inițială Dc exemplu, fie f(x) din figura dc mai jos și aproximația inițială xfU) Aplicând această metodă soluția numerică va oscila între valorile xfl) și xf2) și nu sc va apropia dc soluția exactă x* [01 x Sc observă că dacă sc serie relația dc recurență sub forma J(x(j)) • x(j+1) = J(x(j))x(j) - f(x(j)) este evident că la fiecare iterație se rezolvă un sistem de ecuații liniare a cărui matrice J(x!j>) și termen liber J(x(J))x(J) - f (x(J)) trebuie recalculate Calculul Jacobianului implică determinarea a n2 valori ale derivatelor în x(i) ceea ce nu este deloc simplu Pentru un circuit rezistiv liniar există, însă, o procedură mult mai simplă dc efectuare a iterațiilor Considerăm la început un exemplu: La fiecare iterație rczistorul ncliniar poate fi înlocuit cu un circuit echivalent serie format dintr-o sursa dc tensiune în serie cu un rezistor La prima iterație sursa arc tensiunea electromotoare EfU) și rezistorul are rezistența R(o:i=tg aA Pentru a doua iterație avem 1G! și R!1>=tg a etc Acest circuit sc numește circuitul echivalent discret La fiecare iterație sc calculează parametrii circuitului echivalent discret și sc rezolvă un circuit liniar Această procedură sc poate aplica unui circuit rezistiv ncliniar oricât dc complicat ar fi acesta Circuitul echivalent discret la iterația k+1 al unui rezistor controlat în curent cu caracteristica и=й(і) conține rezistorul cu rezistența R'") = —\,(kîn serie cu sursa de tensiune^ di |г obținută ca în exemplul precedent Pentru un rezistor controlat în tensiune cu caracteristică i=î(u) circuitul echivalent discret are parametrii G(k) = — I (k) și Dacă rezistorul are caracteristici du'u liniare pe porțiuni, parametrii circuitului echivalent discret sc determină ca în figură: 13 Prin utilizarea circuitului echivalent discret aplicarea metodei Ncwton-Raphson sc simplifică considerabil 4 Circuite care funcționează la semnale mici Fie circuitul redresor monoaltcmanță fără filtru din figura (a) unde caracteristica diodei este dată în figura (b) și caracteristica dc transfer în tensiune, este desenată în figura (c) Dacă Ui(t) = sin o t (fig d) atunci din caracteristica de transfer în tensiune rezultă că u2(t) va avea forma din figura (c) în acest caz semnalul dc ieșire arc o forma diferită față dc cel dc intrare datorită neliniaritații caracteristicii de transfer Dacă u1(t) = 1 + sin cot (fig f) atunci semnalul de ieșire este sinusoidal (fig g) ca și cel dc intrare Spunem că un circuit rezistiv neliniar cu o singură sursă dc semnal (cu c sau is variabil în timp) funcționează la semnale mici dacă orice răspuns are forma semnalului de intrare în cazul în care un singur răspuns nu arc aceeași forma cu semnalul dc intrare spunem că circuitul funcționează la semnate mari Pentru un circuit dat funcționarea la semnale mari sau mici depinde dc parametrii semnalului dc intrare Un rezistor ncliniar dipolar funcționează la semnale mici dacă putem aproxima deplasarea 14 punctului dc funcționare pe caracteristică cu deplasarea pe tangenta în punctul static dc funcționare Rezistorul controlat în tensiune din figura dc mai jos cu caracteristica i=g(u) sc comportă din punct dc vedere al semnalului (variațiile dc tensiuneși curent în jurul punctului static de funcționare) ca un rezistor liniar cu conductanța G = ^g(u) du = tg a care se numește p s f • conductanța echivalentă dc semnal mic și depinde dc punctul static dc funcționare Pentru un rezistor controlat în curent cu caracteristică u=f(i) sc definește rezistența echivalentă dc semnal „ df (i) mic R = di Psf • ‘ Fie un tranzistor modelat ca rezistor multipolar cu caracteristicile: uBE = uBi:(i ,u ), i = i (i ,u ) Pentru a determina parametrii circuitului echivalent de B CE C CB CE semnal mic se dezvoltă u și i în serie Taylor în jurul p s f (determinat de u ,i ,i ,u ) Fiecare tensiune și curent are o componentă constantă (corespunzătoare BE0 B0 C0 CE0 p s f și notată cu indicele 0) și o componentă variabilă în timp (dc semnal), adică: u = u BE + u , i BE0 be B = i + i , i = i B0 b C + i , u = u + u C0 c CE CE0 ce Rezultă: u = u BE BE0 du di I B Ps'J'ib du du I CE P •s •f • Uce i C di = ‘ +—I C0 di 1 B di +—q du I CE PsJ-Uce + л w-^b u sau u = r i + au , be b b ce i = в-i +-ce unde c b r c dui dui di „ BE BE в= C 1 rb = dl \P J a = du di B CE B P s f • 1 ~di CE r c Rezultă circuitul echivalent dc semnal mic: 15 Acest circuit poate fi utilizat dacă punctul de funcționare se deplasează într-o zonă din planul fiecărei caracteristici pentru care parametrii rb, a, (3, rc pot fi considerați constanți în aceeași manieră se poate determina un circuit echivalent de semnal mic pentru orice rezistor neliniar Analiza unui circuit care funcționează la semnale mici se poate face pe un circuit echivalent liniar conform algoritmului: 1 Se determină punctul static de funcționare considerând circuitul excitat de sursele de curent continuu (surse independente cu parametrii constanți în timp), sursele de semnal fiind pasivizate 2 Se determină circuitul echivalent de semnal mic asociat p s f determinat la punctul 1 pentru fiecare rezistor neliniar 3 Se face analiza circuitului echivalent de semnal mic obținut prin interconectarea circuitelor echivalente de semnal mic și a rezistoarelor liniare Acest circuit este excitat numai de sursele de semnal, sursele de curent continuu fiind pasivizate Soluția obținută este componenta variabilă în timp a răspunsului circuitului 4 Se verifică dacă condițiile de funcționare la semnal mic sunt satisfăcute pentru toate rezistoarele neliniare Dacă există un singur rezistor neliniar pentru care aceste condiții nu sunt îndeplinite rezultatul analizei de la punctul 3 este incorect 16 II CIRCUITE DINAMICE 1 Existența și unicitatea soluțiilor Existența si unicitatea soluției unui circuit dinamic este legata dc existenta ecuației dc stare in forma normala Daca ecuația x = f (x,t) exista atunci in literatura matematica sc arata ca: daca f este Lipschitzana (pentru orice xt si x, si orice t Ц/(Xj ,t)~ f (x2, r)|| 0 si ||*|| este norma euclidiană) si daca funcția f(O,t) este uniform mărginită, atunci ecuația dc stare arc o soluție unica pentru orice stare inițiala x0 = x( t0 ) Sc construiește arborele normal al circuitului care conține un număr maxim dc condensatoare și un număr minim dc bobine Circuitul dinamic sc reprezintă ca un multiport rezistiv având elementele dinamice din arbore substituite cu surse dc tensiune și elementele dinamice din coarborc substituite cu surse dc curent Existenta formei normale a ecuației dc stare este legata dc : • existenta si unicitatea soluției multiportului rezistiv pentru orice valori ale parametrilor surselor (care inlocuicsc elementele dinamice) conectate la porti, • existenta matricclor A 1 si A deci existenta unor capacități si inductivitati dinamice ncnulc pentru orice valori ale parametrilor dc control ai acestor elemente unde A = diag [ capacitatilc dinamice ale condensatoarelor din arbore, inductivitatilc dinamice ale bobinelor din coarborc ] si A = diag [ capacitatilc dinamice ale condensatoarelor din coarborc, inductivitatilc dinamice ale bobinelor din arbore] Caracterul Lipschitzian al funcției f nu poate fi dedus din caracterul Lipschitzian al ecuațiilor constitutive ale elementelor dc circuit Din acest motiv in teoria circuitelor interesează condiții dc existenta si unicitate exprimate in funcție dc ecuațiile constitutive ale elementelor dc circuit si dc modul dc interconectare a acestor elemente In cazul circuitelor liniare elementele dinamice sunt liniare deci exista Л-1 si Л * 1 Daca multiportul rezistiv arc o soluție si numai una pentru orice valori ale parametrilor surselor 17 independente conectate la porti, atunci exista ecuația de stare in forma normala x = Ax + Bus + B * //s sau x = f (x, t) Se poate arata că în acest caz f (x, t) este Lipschitziana: ||f (x1, t) — f (x2t )|| = | |A( x1 — x2)|| și exista întotdeauna o constanta к astfel incat ||A(x1 — x2)|| 0 i scade iar pentru u 0 soluția exacta (curba exponențiala] soluția numerica / curba exponențiala лп - xn+1 O metodă numerică dc integrare pentru care eroarea totală descrește o dată cu creșterea timpului este o metodă stabilă O metodă care nu arc această proprietate este numeric instabilă, chiar dacă eroarea locală este mică și descrește în timp în exemplul dc mai sus pentru metoda Eulcr explicită în care 23 x1 = (1 - Ah) x0 x2 = (1 - Ah)x1 = (1 - Ah)2x0 xn = (1 - Ah) x este clar că dacă |1 - Ah\ > 1 alunei xn pentru n ~ și deci pentru pași mari de timp 2 metoda este instabilă Pentru ca metoda să fie stabilă trebuie ca |1 - Ah\ x0 1 se vede că, dacă X> 0, pentru n avem xn 0 și deci metodele sunt stabile, indiferent de mărimea h a pasului de timp La o metodă stabilă mărimea pasului de timp este limitată numai de eroarea locală impusă care depinde de problema studiată La metoda Euler explicită £e = ^h2 (e e = k2 h2 pentru metoda Euler implicită, 3 2 e e = k^h pentru metoda trapezului), pasul fiind limitat și de - se poate întâmpla ca A max lucrând cu valori foarte mici ale lui h numărul mare de pași necesari pentru determinarea soluției să ducă la un timp de calcul nejustificat de mare în cazul unei metode explicite, după alegerea pasului de timp, se pornește de la condiția inițială x(t0) și se calculează x(10 + h), x(t0 + 2h), acoperind tot intervalul de timp care prezintă interes în cazul unei metode implicite la fiecare pas se fac mai multe iterații La prima iterație se folosește o metodă explicită și se obține predictorul xk+1 (0) = xk + hf (xk) Valoarea obținută pentru predictor se introduce în membrul drept al ecuației metodei implicite obținându-se o nouă valoare a lui xk+1(1) (în membrul stâng) care la iterația următoare se introduce din nou în (2) membrul drept și se obține 7 (N-1) (N) s a m d pana ce xk- xk+1 Ldqi cu = —\q si и = ф - bobina controlată în flux i = ftf) T =rdqt cu rdq = fq - condensatorul controlat în tensiune q = f(u) q = cu Cdq = — |q si q = Cdqu - condensatorul controlat în sarcină U = f (q) u = Sdqq cu Sdq = — \q Prin înlocuirea fiecărui element neliniar de circuit cu elementele dinamice de mai sus se obține circuitul echivalent la semnale mici care este un circuit liniar Algoritmul de analiză al unui circuit care funcționează la semnale mici este: 1) Se determină punctul static de funcționare prin analiza circuitului în care toate sursele cu parametri variabili în timp sunt pasivizate, bobinele sunt înlocuite cu rezistențe nule și condensatoarele sunt înlocuite cu rezistențe infinite 2) Se calculează parametrii circuitului echivalent de semnal mic (Rdq, Gdq, Ldq, rdq, Cdq, Sdq ,adq,вdq) pentru fiecare element neliniar 3) Se face analiza circuitului echivalent de semnal mic care conține numai sursele independente variabile în timp (sursele de curent continuu fiind pasivizate) 4) Se face verificarea rezultatelor testând în ce măsură aproximarea de semnal mic este corectă pentru fiecare element neliniar de circuit (dacă “excursia” punctului de funcționare pe caracteristica neliniară poate fi aproximată cu “excursia” pe tangenta în p s f ) 7 Metoda potențialelor nodurilor pentru circuite de curent alternativ Asa cum s-a arătat în paragraful 2 , se preferă comanda în tensiune deoarece mărimea de comandă poate fi scrisă ca o diferență de potențiale Ca urmare, circuitul echivalent cu surse de tensiune comandate în curent al bobinelor cuplate nu este potrivit pentru scrierea ecuațiilor metodei nodale Pentru aceste bobine se poate construi un circuit echivalent cu surse de curent comandate în tensiune în acest scop se rezolvă ecuațiile de funcționare ale bobinelor cuplate: 9 U2=j^^2±jwMIi jn rapOrt cu necunoscutele Ij și I2 30 Exemplu Fie bobinele cuplate cu ecuațiile U1 = 2jli + jl2 si U2 = jlț + 2jl2 dc funcționare Rezolvând acest sistem dc ecuații în raport cu It șij 2 rezultă: U1 j U2 j Ii ==1 + U2 Si I2 ==2 + U1 -1 3 -2 3 3-1 j A j- 2 2 adică ecuațiile următorului circuit echivalent: Aplicație Să sc serie ecuațiile potențialelor nodurilor pentru circuitul 0 5F 0 5H ( e=3i-3; e(t) =42 cos2t i (t) = 2sin(2t + -) s 4 Deoarece cuplajul nu sc poate sparge construim circuitul echivalent cu surse dc curent comandate în tensiune al celor două bobine cuplate Cu notațiile din figură ecuațiile dc funcționare ale acestor două bobine sunt: U1 = 2jI1 + Л2 ,U2 = 2jl2 + jI1 ^cc' circuitul echivalent este cel prezentat în exemplul precedent Schema echivalentă în complex este: Sc observă că avem un circuit rezonant RLC serie cu = j-j+3 =3 și două surse ideale dc tensiune care nu sc pot transforma în surse dc curent Alegem ca potențial dc referință o bornă a uncia dintre aceste surse iar pentru cealaltă introducem necunoscuta suplimentară Rezultă ecuațiile: 31 V5 = 0 У1 = j - У1| - V31 = I4 + ju2 3j 3 3 U2 =- У3 У 2 -У 4 = 3І3 У1 - У3 1 - j 1 1 2 У3 +- + — 1 - j 3 - y2 = 1+j+2u1 3 3 1 - j U1 = У1 - У2 y^2)=-1 - j -14 adică un sistem de 9 ecuații cu necunoscutele У1, , У 5, U1, U2,13,I4 Deci algoritmul de scriere a ecuațiilor potențialelor nodurilor este: • se fac toate transformările posibile ale surselor de tensiune în surse de curent și ale comenzilor în curent în comenzi în tensiune, • se alege potențialul de referință astfel încât cât mai multe potențiale ale nodurilor să poată fi exprimate ca sume de tensiuni electromotoare, • considerând și necunoscutele suplimentare (curenții unor surse de tensiune conectate între alte noduri decât cclc dc la punctul precedent și curenții de commandă) sc scrie sistemul de ecuații: Vj £ Yk ViYk =Z Isk k G j k Gi,j k G j și ecuațiile suplimentare 32 CAPITOLUL 2 UTILIZAREA PROGRAMULUI PSPICE 1 Introducere SPICE (Simulation Program with Integrated Circuit Emphasis) este un program destinat simulării circuitelor electrice PSPICE reprezintă varianta SPICE pentru calculatoare personale și a fost dezvoltat dc MicroSim începând cu 1984 Alte versiuni ale programului SPICE sunt : HSP1CE produs dc Mcta-Softwarc, 1-SP1CE produs dc NCSS Timesharing si 1G-SP1CE produs dc A B Associates , SPECTRE produs dc Cadcncc Design, SPICE-Plus produs dc Valid Logic 1 1 Fereastra de comandă Fereastra dc comandă sc poate deschide din Start-Programs sau cu click pe un fișier dc intrare (dc tip * cir sau PSpicc Circuit File) sau pe un fișier dc rezultate (dc tip * dat sau Pspicc Simulation Rcsults) Programul sc mai poate lansa și din fereastra command cu PSPICE numeJișier de intrare sau dintr-n program scris în C, FORTRAN sau alt limbaj dc programare folosind instrucțiuni dc tip system în fereastra dc comandă există următoarele menu-uri: File - ca în toate aplicațiile WINDOWS, poate deschide un fișier, îl poate tipări, memorează ultimele fișiere accesate și ultimele simulări făcute, poate salva sau poate închide un fișier și poate închide fereastra dc comandă Edit - permite modificările uzuale în fișierul deschis în fereastra dc comandă View - permite vizualizarea fișierelor dc intrare (* cir) și dc ieșire (* out), a ferestrei grafice în care sc prezintă rezultatele simulării, a ferestrei în care sc caracterizează global simularea (dacă sunt erori în fișierul dc intrare, dacă un anume tip dc analiză nu converge, dacă simularea s-a încheiat) Simulation -conține comenzile dc începere a simulării, pauză și oprire a simulării Trace - conține comenzile dc adăugare și ștergere a unei reprezentări grafice și comenzi privitoare la cursor Plot- conține comenzi privitoare la formatul dc afișare cum sunt: modificarea scării pe axele X si Y, modificarea intervalului dc afișare, schimbarea mărimii asociate unei axe, introducerea și ștergerea unei axe Y suplimentare Tools- permite introducerea sau eliminarea unor bare cu icoane care reprezintă diverse comenzi Window-un meniu clasic care permite lucrul cu ferestrele din PSPICE Help- un meniu clasic care permite consultarea documentației și accesul pe internet la site-urilc producătorului 33 1 2 Fișierele de intrare și de ieșire Fișierul de intrare Fișierul dc intrare trebuie să aibă extensia CIR Arc o structură liberă, în sensul că nu există o succesiune strictă a liniilor, cu excepția primei linii care reprezintă titlul circuitului si arc valoarea unui comentariu (primul caracter este *) și a ultimei linii ( END) care indică sfârșitul fișierului dc intrare Continuarea unei linii sc face tastând + în prima coloană a liniei următoare PSPICE nu face diferența între litere mari și litere mici Un nume trebuie să înceapă cu o literă și conține maxim 8 caractere Un număr poate fi întreg sau un număr în virgulă mobilă și poate fi însoțit dc un exponent întreg Există următorii factori dc scală: T=1E12 ; G=1E9 ; MEG=1E6 ; K=1E3 ; M=lE-3 ; U=lE-6 ; N=lE-9 ; P=1E-12 ; F=1E-15 Literele care sc succed imediat după un număr și nu sunt factori dc scală sunt ignorate Descrierea circuitului (topologia și descrierea elementelor dc circuit) sc face cu ajutorul liniilor elementelor dc circuit în legătură cu acestea trebuie făcute câteva precizări și anume: * sc utilizează regula dc la receptoare pentru sensurile dc referință ale tensiunilor și curcnților u * nodul 0 reprezintă masa și este obligatoriu, celelalte noduri fiind numerotate succesiv, în ordine crescătoare, începând cu 1 Fiecare nod trebuie să aibă o calc către masă în regim dc curent continuu deoarece în general, chiar dacă sc cerc o analiză în curent alternativ sau în regim tranzitoriu, PSPICE face la început un calcul al punctului static dc funcționare * circuitul nu trebuie să conțină bucle (secțiuni) formate exclusiv din generatoare ideale dc tensiune (curent) Analiza circuitului (tipul analizei și diverșii parametri) este descrisă dc liniile dc comandă Acestea au un punct în prima coloană Fișierul de ieșire Fișierul dc ieșire arc extensia OUT și este generat automat Acesta conține diverse informații ca : listarea sau nu a fișierului dc intrare, listarea sau nu a modelelor, a diverșilor parametri (în funcție dc opțiunile din linia dc comandă OPTIONS), listarea valorilor variabilelor dc ieșire (indicate în linia dc comandă PRINȚ), ctc 1 3 Elementele de circuit PSPICE simulează circuite electrice care conțin următoarele elemente dc circuit, liniare sau ncliniarc: - rezistor R*** - condensator C*** - bobină l*** - bobină cuplată K*** - linie dc transmisie T*** - surse independente 34 - dc tensiune V*** - dc curent I*** Formele dc undă în regim tranzitoriu pot fi dc tipul: - sinusoidal SIN - exponențial EXP - impuls PULSE - liniar pe porțiuni PWL - sinusoidal modulat în frecvență SFFM - surse comandate - sursa dc tensiune comandată în tensiune E*** - sursa dc tensiune comandată în curent II*** - sursa dc curent comandată în tensiune G*** - sursa dc curent comandată în curent F*** Comanda poate fi dc tipul: - liniară - polinomială POLY - o funcție analitică VALUE - dată prin puncte TABLE - o expresie în domeniul Laplacc LAPLACE - dispozitivele semiconductoare uzuale acceptate sunt: - dioda D*** - tranzistorul bipolar q*** - tranzistorul cu efect dc câmp JFET J*** -tranzistorul MOSFET M*** - tranzistorul GaAsFct B*** - subcircuitul x*** (care este definit cu ajutorul elementelor dc circuit prezentate mai sus) Descrierea în ordine alfabetică a elementelor dc circuit este dată în paragraful 2, iar a modelelor asociate lor în paragraful 3, la linia dc comandă MODEL 1 4 Tipurile de analiză Analiza în curent continuu DC - calculează punctul static dc funcționare al circuitului (în curent continuu) OP - determină punctul static dc funcționare Această analiză este automat efectuată înainte dc analiza în regim tranzitoriu, dacă nu sunt specificate condițiile inițiale, precum și înainte dc analiza în regim sinusoidal, pentru determinarea modelelor dc semnal mic ale dispozitivelor ncliniarc dc circuit TF - determină orice funcție dc circuit în curent continuu considerând modelul dc semnal mic corespunzător punctului static dc funcționare 35 SENS - determină scnzitivitățilc unor tensiuni și/sau curcnți în funcție dc toți parametrii modelului dc semnal mic corespunzător punctului static dc funcționare Analiza în regim sinusoidal de semnal mic AC- analizează circuitul în regim sinusoidal, într-o bandă dată dc frecvențe, folosind modelul liniar dc semnal mic corespunzător punctului static dc funcționare NOISE - calculează zgomotul produs la o anumită icșrc (tensiune sau curent) într-o bandă dc frecvențe specificată în comanda asociată AC Analiza în regim tranzitoriu TRAN - integrează ecuațiile modelului dc semnal marc al circuitului pe un interval dc timp dat Condițiile inițiale sunt calculate printr-o analiza în curent continuu în care sc consideră toate sursele independente cu valorile dc la momentul inițial sau sunt date dc utilizator FOUR -face analiza armonică a unei tensiuni sau a unui curent obținut cu comanda TRAN pentru o frecvență fundamentală dată și un număr specificat dc componente armonice Aceste comenzi sunt descrise în detaliu în paragraful 3 1 5 Ecuațiile circuitului în curent continuu și în curent alternativ sc utilizează ecuațiile metodei potențialelor nodurilor Analiza în regim tranzitoriu folosește modelele companion ale elementelor dinamice corespunzătoare unei anumite metode numerice; pentru circuitul rezistiv care rezultă sc utilizează ecuațiile metodei potențialelor nodurilor 1 6 Metodele numerice A Metode de integrare numerică Metodele dc integrare numerică folosite dc PSPICE sunt: - metoda trapezului, - metoda Gear dc ordin 2 până la 6 Unele variante dc PSPICE permit utlizatorului să aleagă între metoda trapezului și metoda Gear dc un anumit ordin maxim La alte variante metoda dc integrare numerică este aleasă dc program B Metode de rezolvare a circuitelor rezistive neliniare' PSPICE utilizează în acest scop metoda Newton-Raphson atât pentru analiza în curent continuu cât și pentru analiza circuitelor care conțin modele companion 1 7 Criterii de convergență a iterațiilor Newton-Raphson Iterațiile Newton-Raphson sc opresc dacă sunt îndeplinite toate condițiile următoare: 36 -eroarea absolută între valorile oricărui curent la ultimele două iterații este mai mică decât valoarea parametrului ABSTOL(dacă nu sc specifică altfel ABSTOL=lpA) -eroarea absolută între valorile oricărei tensiuni la ultimele două iterații este mai mică decât valoarea parametrului VNTOL(dacă nu sc specifică altfel VNTOL=1 pV) -eroarea relativă între valorile oricărui curent sau oricărei tensiuni la ultimele două iterații este mai mică decât valoarea parametrului RELTOL(dacă nu sc specifică altfel RELTOL =0 1%) Atât valorile parametrilor ABSTOL, VNTOL, RELTOL cât și numărul maxim dc iterații Ncwton-Raphson pot fi schimbate în linia dc comandă OPTIONS Dacă procedeul iterativ nu este convergent sc tipăresc potențialele nodurilor pentru ultima iterație 1 8 Surse de erori Existența unor erori în descrierea circuitului, ca și a unor erori dc calcul este semnalată în fereastra în care sc caracterizează global simularea Descrierea detaliată a acestor erori sc găsește în fișierul dc ieșire Dacă circuitul arc un nod pentru care nu există o calc dc rezistență finită spre masă în curent continuu punctul static dc funcționare nu sc poate calcula și programul sc oprește semnalând această eroare Dacă cele două ieșiri ale liniei dc transmisie nu au între ele o calc dc rezistență finită în curent continuu sau dacă nodurile dc comandă ale unei surse comandate în tensiune nu au între ele o calc dc rezistență finită în curent continuu efectul este același Acest tip dc erori sc elimină conectând corespunzător un rezistor cu rezistența foarte marc Pot exista erori datorate existenței unor bucle formate numai din surse dc tensiune sau unor secțiuni formate numai din surse dc curent Aceste erori sc corectează introducând în bucla respectivă un rezistor cu rezistența foarte mică sau introducând în secțiunea respectivă un rezistor cu rezistența foarte marc Pot apărea probleme dc convergență a iterațiilor Ncwton-Raphson în curent continuu sau în regim tranzitoriu în curent continuu sc recomandă utilizarea comenzii NODESET pentru alegerea unei alte aproximații inițiale Tot în acest caz sc poate încerca utilizarea opțiunii STEPGMIN care generează un proces iterativ dc corectare a conductațci minime din circuit; această opțiune poate fi activată prin comanda OPTIONS în regim tranzitoriu cauza poate fi utilizarea unor modele ncrcalistc, fără capacități și rezistențe parazite (dc exemplu circuite conținând diode și bobine, fără capacități și rezistențe parazite), sau setarea nccorcspunzătoarc a diferitelor erori utilizate drept criterii dc convergență ale iterațiilor Ncwton-Raphson 37 2 Descrierea elementelor de circuit Pspicc utilizează următoarele clemente dc circuit, în ordine alfabetică: B*** - Tranzistor GaAsFET Forma generală: B [ ] , , - nodurile corespunzătoare drcnci, grilei și sursei = numele modelului = factor dc arie Exemple: BIN 100 1 0 GFAST В13 22 14 23 GNOM 2 0 C*** - Condensator Forma generală: C [ ] [IC= ] , = nodurile pozitiv și negativ = valoarea capacității în F IC= = tensiunea inițiala, în V, pe condensator (în acest caz sc folosește comanda TRAN cu opțiunea SKIPBP) Exemple: CLOAD 15 0 20pF CFDBK 3 33 CMOD 10pFIC=1 5v D*** - Diodă ^tnodț^ ^-nod^ Forma generală: D [ ] , = nodurile pozitiv și negativ = numele modelului = factor dc arie Exemple: DCLAMP 14 0 DMOD D13 15 17 SWITCH 1 5 E*** -Sursă de tensiune comandată în tensiune Uc 38 Forma generală: E E POLY( ) >* + >* E VALUE={ } E TABLE { } * E LAPLACE { } { } E FREQ { } * , = nodurile pozitiv și negativ ale porții comandate , = nodurile pozitiv și negativ ale porții dc comandă = u ) > * - sursa este comandată polinomial = dimensiunea polinomului (1 ;2;3) = coeficienții polinomului Exemple: - polinom dc dimensiune 1 c = a0 + afu + a2uc2 + - polinom dc dimensiune 2 c = a0 + ajUcl + a2uc2 + a3ucl2 + a4ucluc2+ a5uc22 + a6ucl3 + + a7ucl2uc2+ a8ucluc22 + a9uc23 + - polinom dc dimensiune 3 c = a0 + ajUcl + a2uc2 + a3uc3+ a4ucl2 + a5ucluc2+ a6ucluc3 + + a7uc22 + a8uc2uc3+ a9uc32+a10ucl3 + anucl42+ + +a13ucluc22 + a^u^u^u^H- 15ucluc32 + + ai6Uc2 + a17uc2 uc3 + a18uc2uc3 +a19uc3 + VALUE= { } - sursa este o expresie analitică de variabilele xk (tensiuni, curcnți prin surse ideale dc tensiune, timp) c=f(xk), k=l,2, Funcții utilizabile în expresii: Funcția Semnificația ABS(x) |x| SQRT(x) 1/2 x EXP(x) ex LOG(x) ln x LOG10(x) log x PWR(x,y) |x|y PWRS(x,y) +|x|y (x>0), -|x|y (x } * - mărimea dc comandă este exp - o expresie care depinde dc o tensiune din circuit, - dependența este exprimată prin puncte dc coordonatc * aranjate în ordinea crescătoare a valorilor , unde valorile sunt ale lui iar valorile sunt ale tensiunii sursei, - continuitatea acestei dependențe este asigurată prin interpolare liniară, iar în afara intervalului valorilor tensiunea sursei ia valori constante corespunzătoare capetelor intervalului LAPLACE { } { } - mărimea dc comandă este exp - o expresie care depinde dc o tensiune din circuit, - factorul dc transfer este dat ca funcția dc variabila complexă s FREQ { } * - mărimea dc comandă este exp - o expresie care depinde dc o tensiune din circuit, - factorul dc transfer este dat pentru anumite frecvențe prin valoarea modulului (în dB) și fazei (în grade) Valorile pentru frecvențele intermediare sc obțin prin interpolare liniară (pentru fază) sau logaritmică (pentru modul) Exemple: EBUFF 1 2 10 11 1 0 EAMP 13 0 POLY(1) 26 0 500 ENLIN 100 101 POLY(2) 3 0 4 0 0 0 13 6 0 2 0 005 ESQRT 10 0 VALUE = {SQRT(V(5))} ETAB 20 5 TABLE {10*V(2)}(-5v,5v) (0v,0v) (5v,-5v) E1POLE 10 0 LAPLACE {V(1)} {1 / (1 + s)} EATTEN 20 0 FREQ {V(100)} (0,0,0 10,-2,-5 20,-6,-10) F*** - Sursă de curent comandată în curent Forma generală: F F POLY( ) >* >* , = nodurile pozitiv și negativ = sursa dc tensiune al cărei curent reprezintă mărimea dc comandă = is/ic , când sursa este comandată liniar POLY ( ) - sursa este comandată polinomial (vezi explicațiile pentru sursa dc tensiune comandată în tensiune E***) Exemple: FSENSE 1 2 VSENSE 10 0 FAMP 13 0 POLY(l) VIN 500 FNLIN 100 101 POLY(2) VI v2 0 0 0 9 0 2 0 005 40 G*** - Sursă de curent comandată în tensiune Forma generală: G G POLY( ) >* + >* G VALUE={ } G TABLE { }= * G LAPLACE { } } G FREQ { } * Explicațiile sunt similare cu cclc dc la sursa E*** (sursă dc tensiune comandată în tensiune) Exemple: GBUFF 1 2 10 11 1 0 GAMP 13 0 POLY(1) 26 0 500 GNLIN 100 101 POLY(2) 3 0 4 0 0 0 13 6 0 2 0 005 GSQRT 10 0 VALUE = {SQRT(V(5))} GTAB 20 5 TABLE{V(2)} = (-5v,5v) (0v,0v) (5v,-5v) G1POLE 10 0 LAPLACE {V(1)} {1 / (1 + s)} GATTEN 20 0 FREQ {V(100)}(0,0,0 10,-2,-5 20,-6,-10) H*** - Sursă de tensiune comandată în curent Forma generală: H H POLY( ) >* >* Explicațiile sunt similare cu cclc pentru sursa F*** (sursa dc curent comandată în curent) Exemple: HSENSE 1 2 VSENSE 10 0 HAMP 13 0 POLY(l) VIN 500 HNLIN 100 101 POLY(2) VI v2 0 0 0 9 0 2 0 005 I*** - Sursă independentă de curent 41 Forma generală: I [[DC] ] [AC [ ]] + [ ] , = nodurile pozitiv și negativ [DC] ] = semnal de curent continuu de tip i = AC [ ] = semnal de curent alternativ de tipul i = * sin(2 *3 14*f*t), unde f este frecvența specificată în linia dc comandă AC = semnal tranzitoriu Semnalele tranzitorii pot fi dc tipul: - exponențial EXP(ij irk rdelay rtc fdelay ftc) i = i 1 + (ipk-i 1){ l-exp[-(t-rdelay)/rtc]}, rdelay [ ] = nodurile corespunzătoare drenei, grilei si sursei = numele modelului = factor de arie Exemple: JIN 100 1 0 JFAST J13 22 14 23 JNOM + 2 0 K*** - bobine cuplate Forma generala: K L >* K >* [ ] L * = numele bobinelor cuplate M = factor de cuplaj 0 = numele modelului Se poate defini bobina, neliniara, cu miez feromagnetic, cu histerezis = factor de scala pentru aria transversala a bobinei neliniare Exemple: KTUNED L3OUT L4IN 8 KXFR1 LPRIM LSEC 99 KXFR2 L1 L2 L3 L4 98 KPOT 3C8 L*** bobina О->-* • • 1 -o (+nodț^ ^-node^ Forma generala: L [model] [IC= ] , = nodurile pozitiv si negativ 45 [model] = numele modelului = valoare in H IC= = valoare inițială pentru curentul din bobina (se folosește când în linia dc comanda TRAN sc specifica SKIPBP) Exemple: LLOAD 15 0 20mH L2 1 2 2e-6 LCHOKE 3 42 LMOD 03 LSENSE 5 12 2uH IC=2mA M ***- tranzistor MOSFET Forma generala: M [L= ] [W= ] +[AD= ] [AS= ] [PD= ] [PS= ] + [NRD= ] +[NRS= ] [NRG= ] [NRB= ] , , = nodurile corespunzătoare drenei, grilei, sursei și substratului = numele modelului L,W = lungimea și lățimea canalului AD,AS = ariile dc difuzare ale drenei și sursei PD,PS = parametrii joncțiunii NRD,NRS,NRG,NRB = numerele dc pătrate echivalente ale difuziilor drenei, sursei, grilei și substratului Exemple: M1 14 2 13 0 PNOM L=25u W=12u M13 15 3 0 0 PSTRONG M2A 0 2 100 100 PWEAK L=33u w=12u AD=288p AS=288p PD=60u PS=60u + NRD=14 NRS=24 NRG=10 Q*** - tranzistor bipolar Q [ ] [ ] = nodurile corespunzatoare colectorului si bazei = nod al substratului ( optional ) = numele modelului = factor de arie Exemple: 46 Q1 14 2 13 PNPNOM Q13 15 3 0 1 NPNSTRONG 1 5 Q7 VC 5 12 [SUB] LATPNP R*** - rezistor l(R] R O——I I -O V( ] Forma generala: R [ ] , = nodurile pozitiv si negativ = numele modelului = valoarea rezistentei in Q Exemple: RLOAD 15 0 2k R2 1 2 2 4e4 R*** N1 N 2 S*** - intrerupator controlat in tensiune f[uc] Forma generala: S , = nodurile pozitiv și negativ ale Întrerupătorului , = nodurile pozitiv și negativ ale tensiunii dc commandă Uc = numele modelului Exemple: S12 13 17 2 0 SMOD SRESET 5 0 15 3 RELAY T*** - linie de transmisie o— - A-± ^ Forma generală: T Z0= [TD= | F= [NL= ]] 47 , = nodurile porților 1 și 2 Z0 = impedanța characteristică TD = întârzierea transmisiei F = frecvența NL = lungimea electrică a liniei, unde lungimea fizică L=NL* Л și Л este lugimea dc undă a semnalului dc frecvență F (valoarea predefinită = 0 25) Exemple: T1 1 2 3 4 Z0=220 TD=115ns T2 1 2 3 4 Z0=50 F=5MEG NL=0 5 V*** - sursa independentă de tensiune , Forma generală: V [[DC] ] [AC [ ]] + [ ] Explicațiile sunt similare cu cele de la sursa independentă de curent I*** Exemple: VBIAS 13 0 2 3mV VAC 2 3 AC 001 VACPHS 2 3 AC 001 90 VPULSE 1 0 + PULSE(-1mV 1mV 2ns 2ns 2ns 50ns 100ns) V3 26 77 DC 002 AC 1 + SIN( 002 002 1 5MEG) W*** - Întrerupător controlat în current Forma generală: W , = nodurile pozitiv și negativ ale întrerupătorului = numele sursei dc tensiune al cărei curent reprezintă mărimea dc commandă = numele modelului Exemple: W12 13 17 VC WMOD WRESET 5 0 VRESET RELAY X*** - subcircuit 48 Forma generala: X [ ]* [PARAMS: = *>] * = nodurile la care este conectat subcircuitul = numele subcircuitului Acesta este definit asfel: SUBCKT [ ]* ENDS Examples: X12 100 101 200 201 DIFFAMP XBUFF 13 15 UNITAMP 3 Liniile de comandă AC - Analiza în curent alternativ Forma generală: AC [LIN][OCT][DEC] [LIN], [OCT], [DEC] - se variază frecvența liniar, pe octave sau pe decade = numărul dc puncte (total pentru varianta LIN, pe octavă pentru varianta OCT și pe decadă pentru varianta DEC) în care sc efectuează analiza AC = frecvența inițială = frecvența finală Example: AC LIN 101 10Hz 200Hz AC OCT 10 1KHz 16KHz AC DEC 20 1MEG 100MEG DC - Analiza în curent continuu Forma generală: DC [LIN] [ ] DC [OCT] [DEC] [ ] DC LIST * [ ] = numele generatorului independent de tensiune sau curent a cărui valoare variază [LIN], [OCT], [DEC] = variație liniară, pe octave sau pe decade a valorii generatorului , = valorile inițiale și finale pentru valoarea generatorului = numărul dc puncte pentru valorile generatorului pentru care sc face analiza DC = pasul dc incrementare pentru valoarea generatorului = alt generator independent al cărui parametru sc variază LIST = analiza DC pentru o listă dc valori ale valorii generatorului Exemple: 49 DC VIN - 25 25 05 DC LIN I2 5mA -2mA 0 1mA DC VCE 0v 10v 5v + IB 0mA 1mA 50uA DC RES RMOD(R) 0 9 1 1 001 DC DEC NPN QFAST(IS) + 1e-18 1e-14 5 DC TEMP LIST 0 20 27 50 80 DC PARAM RS -1 1 0 1 DC SRCNAM VSTART VSTOP VINCR END - Sfârșitul fișierului de intrare PSPICE Forma generală: END ENDS - Sfârșitul zonei de definiție a unui subcircuit Forma generală: ENDS [ ] = numele subcircuitului Exemplu: ENDS 741 FOUR -Analiza Fourier Forma generală: FOUR * = frecvența componentei fundamentale * = variabilele pentru care sc dorește analiza Fourier; această analiză sc poate face numai plecând dc la datele calculate cu analiza în regim tranzitoriu Exemplu: FOUR 10KHz v(5) v(6,7) IC - Conditii initiale pentru regimul tranzitoriu Forma generala: IC = >* = - se atribuie nodurilor valorile inițiale Exemplu: IC V(2)=3 4 V(102)=0 INC - Includerea unui fisier Forma generală: INC = numele fișierului, inclusiv calea Exemple: 50 INC SETUP CIR INC C:\\PSLIB\\VCO CIR LIB - Utilizarea unei biblioteci de modele si subcircuite Forma generală: LIB [ ] = numele librăriei, inclusiv calea Dacă lipsește, sc utilizează biblioteca principală NOM LIB Exemple: LIB LIB OPNOM LIB LIB C:\\PSLIB\\QNOM LIB MC -Analiza Monte Carlo Forma generală: MC [DC][AC][TRAN] * = numar de simulari Monte Carlo [DC] [AC] [TRAN] = tipul de analiză = marimile de ieșire pentru care se face analiza Momte Carlo = funcția aplicată valorilor mărimilor dc ieșire pentru a obține o singură valoare carcatcristică Acestea pot fi : YMAX, MAX, MIN, RISEEDGE(valuc), FALL EDGE(valuc) = diverse opțiuni care se pot seta: [LIST], PLOT, PROBE, ALL, FIRST , EVERY , RUNS Exemple: MC 10 TRAN V(5) YMAX MC 50 DC IC(Q7) MIN LIST MC 20 AC VP(13,5) RISE EDGE(1 0) LIST OUTPUT ALL WCASE -Wort Case Analysis Forma generală: WCASE * Exemple: WCASE DC V(4,5) YMAX WCASE TRAN V(1) FALL EDGE(3 5v) VARY BOTH BY RELTOL DEVICES RL MODEL - Model Forma generala: MODEL [ = [ ]]* = numele modelului = tipul dispozitivului Tip c*** CAP condensator L*** IND bobina R*** RES rezistor D*** D dioda Q*** NPN tr bipolar NPN Q*** PNP tr bipolar PNP 51 Q*** LPNP tr lateral PNP J*** NJF JFET canal N J*** PJF JFET canal P M*** NMOS MOSFET canal N M*** PMOS MOSFET canal P B*** GASFET GaAsFet canal N K*** CORE miez neliniar s*** VSWITCH întrerupător controlat în tensiune W*** ISWITCH întrerupător controlat în curent N*** DINPUT dispozitiv digital de intrare O*** DOUTPUT dispozitiv digital de iesire = [ ] - setarea anumitor parametri la valorile și eventual cu toleranțele în procente Urmează lista parametrilor câtorva tipuri dc dispozitive și anume: ** pentru rezistor R*** Nume Parametru Unitate Val predefinita R coeficient dc amplificare pentru rezistența 1 0 TC 1 coeficient liniar dc temperatura 0C-1 0 0 TC 2 coeficient pătratic dc temperatura 0C-2 0 0 TCE coeficient exponențial dc temperatura 0/0 - 0C 0 0 = RNOM* R * [ 1+ TC1( T - TNOM) + TC2( T - TNOM)2] sau = RNOM * R* 1 01TCE ( T - TNOM ) RNOM = valoarea rezistenței la temperatura TNOM TNOM = temperatura nominală Poate fi specificată în linia dc commandă OPTIONS ** pentru condensator C*** Nume Parametru Unitate Val predefinita C coeficient de amplificare pentru capacitate - 1 0 VC1 coeficient liniar de tensiune V-1 0 0 VC2 coeficient patratic de tensiune V-2 0 0 TC1 coeficient liniar de temperatureă 0C-1 0 0 52 TC2 coeficient patratic de temperatureă 0C-2 0 0 = CNOM C (1+VC1 r3 +VC2 0 2) CNOM = valoarea capacității la temperatura TNOM TNOM = valoarea nominală a temperaturii Se poate specifica în linia de comandă OPTIONS 0= tensiunea pe condensator T = temperatura ** pentru bobina L*** Nume Parametru Unitate Val predefinita L coeficient dc amplificare pentru inducțantc - 1 0 IL1 coeficient liniar dc curent A-1 0 0 IL2 coeficient pătratic dc curent A-2 0 0 TC1 coeficient liniar dc temperatura 0C-1 0 0 TC2 coeficient pătratic dc temperatura 0C-2 0 0 = LNOM*L*(1+IL1 i+IL2 i2) LNOM = valoarea inductanței la temperatura TNOM TNOM = temperatura nominală i = curentul prin bobbină T = temperatura ** pentru întrerupător controlat în tensiune S*** Nume Parametrii Unitate Val predefinita RON Rezistența “on” O 1 0 ROFF Rezistența “off” O 1E+6 VON tensiunea dc control pentru starea on V 1 0 VOFF tensiunea dc control pentru starea “off’ V 0 0 * * pentru întrerupător controlat în curent W*** Nume Parametri Unitate Val predefinita RON Rezistența “on” O 1 0 ROFF Rezistența “off’ O 1E+6 ION curent dc control pentru starea “on” A 1E-3 IOFF curent dc control pentru starea “off’ A 0 0 53 ** pentru dioda D*** Nume Parametri Unitate Val predefinită IS curent dc saturație A 1 0E-14 RS Rezistența ohmica Q 0 0 N coeficient dc emisie - 1 0 TT timp dc tranzit s 0 0 CJO capacitatea joncțiunii F 0 0 VJ potențialul joncțiunii V 1 0 M coeficientul dc gradare - 0 5 EG energie dc activare cV 1 11 XTt exponent dc temperatura al curentului dc saturație - 3 0 KF coeficient dc zgomot flickcr - 0 AF exponent dc zgomot flickcr - 1 FC coeficient pentru formula capacității în polarizare directă - 0 5 BV tensiunea dc străpungere V IBV curentul la tensiunea dc străpungere A 1 0E-3 Parametrii modelelor altor dispozitive sunt prezentați în detaliu în bibliografic Exemple: MODEL RMAX RES (R=1 5 TC= 02 TC2= 005) MODEL QDRIV NPN (IS=1e-7 BF=30) MODEL DLOAD D (IS=1e-9 DEV 5% LOT 10%) NODESET - Setarea valorilor potențialelor nodurilor în curent continuu pentru îmbunătățirea convergenței iterațiilor Newton-Raphson Forma generală: NODESET = >* = - nodurilor li se atribuie valorile * Exemplu: NODESET V(2)=3 4 V(3)=-1V NOISE - Analiza zgomotului Forma generală: NOISE [ ] = variabila pentru care se analizeaza zgomotul ; această analiză se poate ultiliza numai împreună cu analiza în curent alternativ AC = numele generatorului independent dc tensiune sau curent care constituie sursa dc zgomot = frecvența pentru care sc tipărește suplimentar contribuția fiecărei surse la zgomotul total al circuitului 54 OP - Punct static de funcționare în curent continuu Forma generală: OP OPTIONS - Opțiuni Forma generală: OPTIONS [ *] [ = *] , = opțiuni care nu cer sau cer atribuirea unor valori Acestea sunt: ACCT - sc tipăresc timpi dc calcul EXPAND - indică elementele din subcircuitc și parametrii acestora LIBRARY - listează liniile din fișiere ale bibliotecilor LIST - listează datele dc intrare NODE NOECHO NOMOD NOPAGE OPTS ABSTOL CHGTOL CPTIME DEFAD DEFAS = DEFL = DEFW = GMIN = ITL1 = Valoarea - tipărește tabelul nodurilor - suprimă listarea - suprimă listarea parametrilor modelului - suprimă paginarea automată - tipărește valori ale diferitelor opțiuni = - schimbă valoarea erorii absolute pentru curent Valoarea predefinită este IpA - - schimbă valoarea erorii absolute pentru sarcină/ flux Valoarca predefinită este 1 0E-14 = -fixează valoarea maximă (în secunde) a CPU time = - schimbă valoarea arici dc difuzare a drenei pentru dispozitive MOS Valoarea predefinită este 0 0 - schimbă valoarea arici dc difuzare a sursei pentru dispozitive MOS Valoarea predefinită este 0 0 - schimbă valoarea lungimii canalului pentru dispozitive MOS Valoarea Predefinită este 100 0 țtm - schimbă valoarea lărgimii canalului pentru dispozitive MOS Valoarea predefinită este 100 0 цт - schimbă valoarea lui GMIN, conductanța minimă acceptată a unei laturi Valoarea predefinită este 1 0 E-12 - schimbă numărul maxim admis dc iterații pentru determinarea p s f ITL2 = ITL4 = ITL5 = tranzitoriu; predefinită este 150 - schimbă numărul maxim admis dc aproximații inițiale pentru analiza în curent continuu ; valoarea predefinită este 20 - schimbă numărul maxim admis dc iterații pentru fiecare punct în cadrul analizei regimului tranzitoriu ; valoarea predefinită este 10 - schimbă numărul maxim total admis dc iterații în cadrul regimului valoarea predefinită este 5000, pentru ITL5=0 sc omite această limită LIMPTS = - schimbă numărul maxim dc puncte pentru un grafic din fișierul * out METHOD = - stabilește metoda dc integrare numerică, cc poate fi Gear sau metoda trapczclor va fi Gear sau trapczoidal Valoarca predefinită este trapczoidal 55 MAXORD= - fixează ordinul maxim pentru metoda dc integrare numerică Gear Poate fi între 2 si 6 Valoarea predefinită cstc2 NUMDGT = - schimbă numărul dc cifre semnificative ale rezultatelor numerice trebuie să fie între 0 și 8 Valoarea predefinită este 4 PIVREL = - raportul maxim admisibil între valoarea pivotului unei matrici și valoarea celui mai marc element dintr-o coloană (valoare predefinită l E-3) PIVTOL= - schimbă valoarea minimă absolutcă admisă pentru un pivot; valoarea predefinită este 1 0 E-13 RELTOL = - valoarea erorii relative admise, atât pentru tensiune, cât și pentru curent (valoarea predefinită 0 001) - temperatura nominală (valoarea prcdcfinită27°C) TNOM = VNTOL= - eroarea absolută pentru tensiune (valoarea predefinită ІцѴ) WIDTH = - lățimea fișierului dc ieșire, în număr dc caractere (valoare predefinită 80) Exemple: OPTIONS NOECHO NOMOD RELTOL= 01 OPTIONS ACCT DEFL=12u DEFW=8u PARAM - Parametru global Forma generală: PARAM = >* ={expr}>* = - parametrului i se atribuie valoarea ={expr} - parametrului i se atribuie expresia {expr} Referirea la valoarea unui parametru global se face utilizand acoladele { } Exemple: PARAM pi=3 14159265 PARAM RSHEET=120, VCC=5V PLOT - Trasarea unui grafic Forma generală: PLOT [DC][AC][NOISE][TRAN] [ [ *] [( , )] ]* [DC],[AC],[NOISE],[TRAN] - tipul de analiza pentru care se cere trasarea graficului = numele variabilei de ieșire pentru care se trasează graficul Numărul maxim al acestor variabile este de 8 ( , ) = punct de origine al axelor pentru graficul respectiv Exemple: PLOT DC V(3) V(2,3) V(R1) I(VIN) PLOT AC VM(2) VP(2) VG(2) PLOT TRAN V(3) V(2,3) (0,5V) ID(M2) I(VCC) (-50mA,50mA) PRINT - Tiparirea rezultatelor Forma generală: PRINT [DC][AC][NOISE][TRAN] [ *] 56 [DC],[AC],[NOISE],[TRAN] - tipul de analiză pentru care se cere tipărirea rezultatelor = variabilele care urmează a fi tipărite Exemple: PRINT DC V(3) V(2,3) V(R1) IB(Q13) PRINT AC VM(2) VP(2) VG(5) II(7) PRINT NOISE INOISE ONOISE DB(INOISE) PROBE -Postprocesorul grafic PROBE Forma generală: PROBE[/CSDF] PROBE[/CSDF] [ *] = variabilele pentru care se dorește o analiză graphică PROBE permite vizualizarea unor expresii dc variabile dc ieșire Variabilele dc ieșire posibile sunt: - pentru regimul dc curent continuu sau tranzitoriu: V( ) = potențialul nodului V( , ) = tensiunea între nodul și V( ) = tensiunea pe clementul dipolar Vx ( ) = potențialul nodului x al unui clement Vxy( ) = tensiunea între nodurile x și у ale unui clement multipolar I( ) = curentul prin clementul dipolar Ix( ) = curentul prin borna x a unui clement multipolar poate fi numele oricărui clement dc tipul: C/D/E/F/G/H/I/L/R/V xy pot fi oricare din bornele: - D/G/S (B) - D/G/S (J) - D/G/S/B (M) - C/B/E/S (Q) - pentru regimul sinusoidal: M - amplitudinea DB - amplitudinea in dB P- faza G - întârzierea de grup R - partea reală I - partea imaginară - pentru analiza zgomotului: INOISE - zgomotul dc intrare ONOISE - zgomotul dc ieșire DB (INOISE) -zgomotul dc intrare în dB DB (ONOISE) - zgomotul dc ieșire în dB Funcțiile utilizabile sunt: 57 Funcție Semnificație ABS(x) (x) SGN(x) +l(x>0),0(x=0),-l(x * * = variabilele dc ieșire pentru care sc dorește analiza senzitivităților în raport cu toți parametrii circuitului Exemplu: SENS V(9) V(4,3) I(VCC) STEP - Analiza în pași 58 Forme generale: STEP [LIN] STEP [OCT][DEC] STEP LIST * [LIN], [OCT],[DEC], LIST - se variaza variabila liniar, pe octave, pe decade sau în punctele * specificate în listă = valorile dc început și sfârșit ale variabilelor = incrementul = numărul de puncte Exemple: STEP VIN - 25 25 05 STEP LIN I2 5mA -2mA 0 1mA STEP RES RMOD(R) 0 9 1 1 001 STEP TEMP LIST 0 20 27 50 80 STEP PARAM X 1 5 0 1 SUBCKT - Definirea unui subcircuit Formele generale: SUBCKT [ *] [PARAMS: [= ]* ] = numele subcircuitului * = nodurile extreme ale subcircuitului Niciunul nu poate fi zero PARAMS: [ = ]* - setarea unor parametri Exemple: SUBCKT OPAMP 1 2 101 102 SUBCKT FILTER IN OUT PARAMS: CENTER, WIDTH=10KHz TEMP - Temperatura Forma generala: TEMP * * = valorile temperaturii pentru care se dorește analiza circuitului Exemple: TEMP 125 TEMP 0 27 125 TF - Calculul unei funcții de transfer în curent continuu Forma generală: TF = variabila dc ieșire = numele generatorului a cărui valoare reprezintă variabila dc intr Exemple: TF V(5) VIN 59 TF I(VDRIV) ICNTRL TRAN - Analiza de regim tranzitoriu Forma generală: TRAN[/OP] [ [ ]] [SKIPBP] = pasul de timp pentru care se afișează rezultatul = timpul final pentru care sc dorește analiza tranzitorie = timpul inițial pentru care sc dorește afișarea rezultatelor = limita superioara a pasului de timp; daca nu se specifica se considera ceiling=ftime/50 pentru un circuit care contine elemente dinamice sau ceiling=pstep pentru un circuit rezistiv [SKIPBP] = folosirea condițiilor inițiale definite in liniile diverselor elemente dc circuit sau calculate folosind setările din linia dc comanda IC Exemple: TRAN 1nS 100nS TRAN/OP 1nS 100nS 20nS SKIPBP TRAN 1nS 100nS 0nS 1nS WATCH - Furnizarea rezultatelor în forma numerică în timpul analizei circuitului Forma generală: WATCH [DC] [AC] [TRAN] >* [DC], [AC], [TRAN] = tipul analizei = numele variabilei (până la trei variabile ) și limitele sale inferioară și superioară Exemple: WATCH DC V(1) (1,5) I (R5) WATCH TRAN VC(2) (-0 8, 0 8) IL(L) (0 1mA, 0 6mA) WIDTH - Lungimea liniilor în fișierul de ieșire Forma generală: WIDTH OUT= = numărul dc coloane fixat pentru liniile din fișierul dc ieșire Poate fi între 80 și 133 Valoarea predefinită este 80 Exemple: WIDTH OUT=80 60 CAP 3 LUCRĂRI DE LABORATOR LUCRAREA 1 CIRCUITE LINIARE DE CURENT CONTINUU 1 PROBLEME Sa sc simuleze cu PSPICE următoarele circuite si să sc explice rezultatele: Ej = 9V; E4=9U3 ; E6 = 12V; Is3 = IA; r, - - 6 O : fi = 3 Q ; r5 = 3 Q ; Fig 1 1 E^V; E2 = 2*Ii; a) E3 = IV; Isl = lA;fi=r2=lQ b) E4 = 2V; Tj = r2 = IQ; r3 = 2Q; c) r4 = 2Q; Fig 1 2 2 INTERPRETAREA REZULTATELOR Pr 1 1 Sc observa ca circuitul nu arc bucle formate numai din surse dc tensiune sau secțiuni formate numai din surse dc curent Circuitul liniar arc deci o soluție unică Soluția determinata cu SPICE sc verifica cu bilanțul puterilor utilizând datele din fișierul * out Pr 1 2 In cazul unui circuit liniar cu surse comandate sc poate ca determinantul matricei sistemului dc ecuații al circuitului sa sc anuleze pentru anumite valori ale parametrilor dc commandă Pentru a determina existenta și unicitatea soluției circuitului din figura 1 2 pentru cazurile a), b) și c) sc determina generatorul echivalent dc tensiune al circuitului dintre bornele A șiB A В Fig 2 1 Rezistenta echivalenta între bornele circuitului pasivizat, RAB0 poate fi determinată după cum urmează Deoarece acest circuit conține o sursa comandata, pentru calculul lui RAB0 se aplica între A si В un curent dc 1 A ca în figura 2 1 Rezulta tensiunea U din care se determină RAB0 = U Sc determina generatorul echivalent dc tensiune intre bornele AB Pentru determinarea rezistentei echivalente R^ sc pasivizează circuitul 61 Rezulta ecuațiile h +1 з = h ■ 21 j — 2I3 =0 213 = 217 + U Cu I2 = IA Din ultimele ecuații rezulta U=0 respectiv: R4B0 = — = 0^ lA Pentru a determina ил4В0 sc scriu ecuațiile circuitului: -21 =21 — U l 4 1 ABO Rezulta 3 ѵ4=1 + Л =>v4=— A It+It+2It=2 2 6 [ил„,=2Ц+2Ц p»» 7+4-U Rezulta generatorul echivalent dc tensiune la bornele AB a) In acest caz иABO ~E3 = 0 si 01 = 0 si circuitul arc o infinitate dc soluții (orice leR)-, r 62 b) In acest caz U AB0 - E4 * 0 rezulta 0 ■ I * 0 si circuitul nu are soluție; c) In acest caz I U AB0 R4 si pentru R4 * 0 circuitul are o soluție unica 3 INSTRUCȚIUNI SPICE Sc vor utiliza următoarele instrucțiuni: a) Descrierea elementelor de circuit - resistor R [ ] R = numele rezistorului , = nodurile intre care se conectează = numele modelului (opțional) = valoarea rezistentei, in Q - surse independente de curent I [[DC] ] [AC [ ]] [ ] , = nodurile pozitiv si negativ [[DC] ]= semnal de curent continuu de tip I = , in A - surse independente de tensiune V [[DC] ] [AC [ ]] [ ] , = nodurile pozitiv si negativ [[DC] ] = semnal de tensiune continua de tip E = , in V - surse comandate liniar - sursa de tensiune comandata in tensiune (VCVS) 63 E , = nodurile pozitiv si negativ ale sursei , = nodurile pozitiv si negativ ale portii de comanda = E/Uc - sursa de tensiune comandata in curent (CCVS) H , = nodurile pozitiv si negativ = sursa de tensiune al carei curent reprezinta marimea de comanda = E/Ic о ■s-conlrol?- о і 9 «-nodei* E=fțic) - sursa de curent comandata in tensiune (VCCS) G , = nodurile pozitiv si negativ ale sursei , = nodurile pozitiv si negativ ale portii de comanda = Is/Uc O - sursa de curent comandata in curent (CCCS) F , = nodurile pozitiv si negative ale sursei = sursa de tensiune al carei curent reprezinta marimea de comanda = ЬД Q b) Descrierea tipului de analiza de curent continuu OP - realizeaza determinarea punctului de functionare in current continuu si tipareste valorile potentialelor nodurilor si ale curentilor prin sursele de tensiune in fisierul * out DC [LIN] [ ] DC [OCT] [DEC] [ ] DC LIST ] = numele sursei a carei valoare se variaza; [LIN],[OCT],[DEC] = variatie de tip liniar, pe octave, sau decadica a valorii sursei; LIST = variatie dupa o lista de valori ; , = valorile de inceput si sfsrsit ale parametrului considerat; = increment; 64 = numărul de puncte in intervalul considerat; = cea de-a doua sursa a cărei valoare se variază (optional) ; LUCRAREA 2 CIRCUITE NELINIARE DE CURENT CONTINUU 1 PROBLEME 1 1 ) Sa se simuleze cu PSPICE urmatorul circuit : E1= 1 10 V, r1= 1 kQ, r2= 2 kQ Sa se reprezinte grafic caracteristica de transfer v(2) = f(V(1)) si sa se indice utilitatea acestui circuit 1 2 ) Sa se simuleze cu PSPICE urmatorul circuit: r1= 3 0 Q, E1=10 0 V a) Sa se aprecieze numarul de solutii al circuitului prin metoda dreptei de sarcina b) Sa se rezolve circuitul plecand de la urmatoarele valori pentr V(2) (se foloseste comanda NODESET) : 0 0V ; 3 0V ; 7 0V ; 65 2 INTERPRETAREA REZULTATELOR La problema 1 1 se obtine cu SPICE urmatoarea caracteristica de transfer: Se constata ca pentru valori mai mari sau egale cu 6V ale tensiunii de intrare V(1), tensiunea de iesire V(2) este egala cu 4V Rezulta ca circuitul este un stabilizator de tensiune La punctul a) al problemei 1 2 se traseaza dreapta de sarcina cu ecuatia E = Ri + u respectiv: 10 = 3i + u Dreapta se intersecteaza cu axele in punctele A de coordonate (0, 3,33) si B de coordonate (10, 0) Se constata ca dreapta de sarcina se intersecteaza cu caracteristica rezistorului neliniar in punctele corespunzatoare valorilor U1 = 1,1765 V, U2 = 3,2174 V si U3 = 4,6102 V Rezulta ca circuitul are trei solutii, deci un numar impar de solutii in conformitate cu teorema de existenta de la paragraful 2 4 1 2 din curs Programul SPICE rezolva circuitul utilizand metoda Newton-Raphson In cazul in care circuitul are mai multe solutii, prin aceasta metoda se determina numai una dintre ele si anume aceea care este “mai apropiata” de aproximatia initiala In problema, aproximatia initiala se introduce cu instructiunea NODESET Atribuind acestuia valorile 0V, 3V si 7V se obtin succesiv cele trei solutii Totodata se verifica faptul ca daca aproximatia initiala nu este specificata, programul SPICE ii atribuie acesteia valoarea nula V(0) = 0 si in acest caz se determina prima solutie 3 INSTRUCTIUNI SPICE Se utilizeaza urmatoarele linii pentru descrierea circuitului: 66 1) Descrierea unei diode o -t'-3 ° D [ ] , = nodurile pozitiv si negativ = numele modelului = numarul de diode montate in paralel (daca acest parametru nu este precizat se considera area =1); 2) Descrierea unui model MODEL [ = ] = numele modelului ( de exemplu cel declarat in linia de descriere a diodei); = tipul modelului (pentru dioda este D); = - reprezinta setarea unor parametrii ai modelului Principalii parametri ai diodei sunt: Parametru IS RS CJ0 BV IBV Nume curent de saturatie rezistenta serie capacitaea jonctiunii la tensiune nula tensiunea inversa de strapungere curentul la tensiunea BV Unitate de masura A Q F V A Valoare predefinita 1 0E-14 0 0 1 0E-3 3) Descrierea surselor comandate neliniar: - comanda de tip polinom: E POLY( ) >* + >* F POLY( ) >* >* G POLY( ) >* + >* H POLY( ) >* >* = dimensiunea polinomului (1;2;3) = coeficientii polinomului Exemplu: - polinom de dimensiune 1 e = a0 + a1uc + a2uc2 + - comanda de tip o expresie analitica: E VALUE=[ ] 67 G VALUE=[ ] - o expresie analitica; - comanda de tip tabel (intre puncte se considera extrapolarea liniara): E TABLE [ ] * G TABLE [ ]= * Cu acest ultim tip de surse comandate se modeleaza rezistoarele neliniare cu caracteristici liniare pe portiuni (PWL) LUCRAREA 3 CIRCUITE CU EXCITATIE SINUSOIDALA CARE FUNCTIONEAZA LA SEMNALE MICI SAU LA SEMNALE MARI 1 PROBLEME 1 1 ) Sa se simuleze cu PSPICE urmatorul circuit: Rb = 1MQ ; Rc = 4kQ; E1 = 5V; Vcc = 9V; e2 = E2 sin (mt); TB: в= 200 ; Cje = 2pF; Cjc = 2pF; Vje = 0 6V; Vjc = 0 6V; a) Sa se reprezinte grafic caracteristica de transfer in curent continuu Ve = f (Vi), unde Vi e , iar Ve este tensiunea pe rezistorul Rc b) Sa se faca analiza de curent alternativ la semnal mic a circuitului pentru E2=1 0 V si respectiv E2= 20 0 V Sa se reprezinte grafic Ve=f (frecv) , frecventa variind decadic intre 0 1 Hz si 100Mhz Sa se discute valabilitatea rezultatelor obtinute prin analiza de semnal mic pentru cele doua valori ale lui E2, timand seama de caracteristica determinate la punctul a) 2 INTERPRETAREA REZULTATELOR La punctul a) al problemei, din analiza caracteristicii de transfer in curent continuu Ve = f (Vj) prezentata in figura 1, rezulta ca aceasta prezinta doua portiuni liniare distincte: 68 1 0 V v (5 )-v (4 ) V 1 Fig 1 prima porțiune pentru Vi 11,25V Acest rezultat se datoreaza faptului ca tranzistorul are caracteristici neliniare de intrare si de iesire atunci cand functioneaza la semnal mare Prima portiune liniara a carcateristicii ue = f (u,) corespunde conditiilor de functionare la semnal mic La punctul b) al problemei in urma analizei in curent alternativ a circuitului rezulta Ve = f (frecv) pentru E2 = 1 0 V, fig 2 si E2 = 20 0 V, fig 3 V (5 )-V (4 ) F r equ ency Fig 2 69 2 0 V Frequency V (5 )-V (4) Fig 3 Se constata ca alura celor doua curbe este identica, curba din figura 3 reproducând la o alta scara curba din figura 2 Luand in considerare caracteristica determinata la punctul a) rezulta ca rezultatul obtinut la punctul b) pentru E2 = 20 0 V nu este valabil, deoarece pentru Vt > 11,25 V circuitul functioneaza la semnal mare, deci modelul liniar la semnal mic nu este valabil 3 INSTRUCȚIUNI SPICE Se vor utiliza urmatoarele linii: a) descrierea tranzistorului bipolar: Q [ ] = nodul corespunzator colectorului = nodul corespunzator bazei = nodul corespunzator emitorului = nodul corespunzator substratului (optional) = numele modelului b) descrierea modelului tranzistorului bipolar : MODEL [ = ] = numele modelului (declarat in linia de descriere a tranzistorului) =tipul modelului (pentru tranzistor bipolar este PNP sau NPN): = - reprezinta setarea unor parametrii ai modelului Principalii parametri ai tranzistorului bipolar (modelul Gummel-Poon) sunt: Parametru IS BF RB RE RC Nume curent de amplificarea rezistenta rezistenta rezistenta 70 saturatie in curent в bazei emitorului colectorului Unitate de masura A - Q Q Q Valoare predefinita 1 0E-16 100 0 0 0 Parametru Vje Vjc Cje Cjc Nume tensiune de deschidere BE tensiune de deschidere BC capacitatea BE (jonctiune nepolarizata) capacitatea BC (jonctiune nepolarizata) Unitate de masura V V F F Valoare predefinita 0 75 0 75 0 0 c) descrierea analizei in curent alternativ a circuitului echivalent de semnal mic: AC [LIN][OCT][DEC] [LIN][OCT][DEC] = variatia liniara, pe octave, sau decadica a frecventei; = numarul de puncte in intervalul considerat; , = capetele intervalului de frecventa considerat; Pentru sursele independente forma generala este: V AC [ ] I AC [ ] , = nodurile pozitiv si negativ AC [ ] = semnal de curent alternativ de tipul x = * sin (2 nf t + phase), unde x este tensiune sau curent si f este frecventa specificata in linia de comanda AC 71 LUCRAREA 4 CIRCUITE LINIARE DE CURENT ALTERNATIV 1 PROBLEME 1 1 ) Sa se simuleze cu PSPICE urmatorul circuit: Se dau: r1=10Q; r2=1Q; r3=2Q; c2=1 5mF; L1=50mH; L2=25mH; L3=50mH; M=25mH; e(t)=30 sin(2nft + n/2); is(t)=2 sin(2nft); f e (10Hz,200Hz) Se cere sa se calculeze curentul prin x2 (modulul, argumentul, partea reala si partea imaginara) pentru f=50Hz si sa se verifice concordanta acestor rezultate 1 2 ) Sa se traseze cu PSPICE caracteristicile de amplitudine si de faza ale amplificarii in tensiune Au= Ve/e pentru urmatoarele filtre in banda de frecvente 1Hz-1MHz Sa se indice tipul fiecarui filtru r1 =1kQ; c1 = 8nF; c2 =11 5nF L1 = 3mH;L2 =12mH; e(t)=1 sin (mt) 72 Fig b) Гі =1,5 kQ; c1 = 10nF; c2 = 8uF Li = 2mH;L2 =12mH; e(t)=1 sin (wt) Fig c) r1 =1,5 kQ; c1 = 10nF; c2 = 8uF L1 = 2mH; c3=1nF e(t)=1 sin (wt) ; Rs=1 5kQ Fig d) r1 =1,5 kQ; c1 = 0,1uF; c2 = 40uF L1 = 10mH; L2 =20mH; c3=1uF e(t)=1 sin (wt) ; L3=4mH 2 INTERPRETAREA REZULTATELOR La problema 1 1 chestiunile care se vor urmari sunt: - modul de descriere a elementelor de circuit, bobine, bobine cuplate si condensatoare, in circuitele de curent alternativ; - modul de obtinere si tiparire a marimilor electrice, curenti si tensiuni, la o frecventa precizata; La problema 1 2 , pentru fiecare punct, se va trasa caracteristica de amplificare si se vor obtine urmatoarele tipuri de filtre: trece jos, trece sus, trece banda si opreste banda 3 INSTRUCTIUNI SPICE Se vor utiliza urmatoarele linii: a) descrierea bobinelor: L [ ] [IC= ] = nodurile intre care se conecteaza model = numele modelului (optional) 73 value = valoarea inductivitatii (in H) IC= = condiția inițiala (pentru curent) (se folosește pentru TRAN cu opțiunea SKIPBP) b) descrierea cuplajului K L >* L >* = bobinele intre care exista cuplaj = factorul de cuplaj (K [ ] [IC= ] , = nodurile pozitiv si negativ = valoarea capacitatii in F IC= = tensiunea initiala, in V, pe condensator (se foloseste pentru TRAN cu optiunea SKIPBP) d) descrierea analizei in curent alternativ a circuitului echivalent de semnal mic: AC [LIN][OCT][DEC] [LIN][OCT][DEC] = variatia liniara, pe octave, sau decadica a frecventei; = numarul de puncte in intervalul considerat; , = capetele intervalului de frecventa considerat; Pentru sursele independente forma generala este: V AC [ ] I AC [ ] , = nodurile pozitiv si negativ AC [ ] = semnal de curent alternativ de tipul x = * sin (2 nf t + phase), unde x este tensiune sau curent si f este frecventa specificata in linia de comanda AC d) Comanda PRINT pentru analiza de tip AC PRINT AC VM(1) VP(5) IR(r1) VI(5,4) VM(1) = amplitudinea potentialului nodului 1 (modulul numarului complex); VP(5) = faza potentialului nodului 5 (argumentul numarului complex); IR(r1) = partea reala a curentului prin r1; VI(5,4) = partea imaginara a tensiunii intre nodurile 5 si 4 74 LUCRAREA 5 SUBCIRCUITE ANALIZA IN REGIM TRANZITORIU 1 PROBLEME 1 ) Sa se determine caracteristica amplificării A= U2 / U1 pentru urmatorul circuit cu amplificatoare operationale: a) E1=10-3sin(2nft) (V), fe (0 1Hz , 10MHz) Pentru amplificatorul operational se va utiliza urmatorul model, care va constitui un subcircuit: ri = 1MQ; re=0 1mQ; Ce=0 1pF b) sa se simuleze in regim tranzitoriu acelasi circuit cu f=500Hz si E1 = 1mV, E1 = 10V; sa se vizualizeze tensiunea de iesire sis a se explice rezultatele obtinute 2 ) Sa se simuleze in regim tranzitoriu urmatorul circuit : 75 Rb = 1MQ ; R = 4kQ; Ei = 5V; Vcc = 9V; e2 = E2 sin (rot); f=10kHz; TB: в= 200 ; Cje = 2pF; Cjc = 2pF; Vje = 0 6V; Vjc = 0 6V; Sa se faca analiza pentru E2=1 0 V si respectiv E2= 20 0 V Sa se reprezinte grafic Ve Sa se compare formele de unda obtinute pentru cele doua valori ale lui E2 Sa se discute valabilitatea rezultatelor obtinute la lucrarea L3 pentru acelasi circuit 2 INTERPRETAREA REZULTATELOR La problema 1 chestiunile care se vor urmari sunt: - Modul de utilizare a subcircuitelor; - Analiza in curent alternativ ( AC) excitand circuitul in gama de frecvente 0,1 Hz - 10 Mhz; - Se va verifica prin analiza in regim tranzitoriu ca pentru amplitudinea de 1 mV la intrare se obtine acelasi rezultat ca la analiza in curent alternativ deci circuitul functioneaza la semnale mici Daca tensiunea de intrare are amplitudinea de 10V forma de unda la iesire nu este sinusoidala deci circuitul functioneaza la semnal mare si rezultatele obtinute cu analiza de semnal mic (AC) nu sunt valabile La problema 2 simuland circuitul in regim tranzitoriu se obtin pentru E2=1 0 V figura 1 si pentru E2= 20 0 V figura 2 V(5)-V(4) * V(2) Time Fig 1 76 4 0 V V(5)-V(4) V(2) X o Time Fig 2 Din figura 2 rezulta ca forma de unda a tensiunii Ve pe rezistorul de sarcina Re, care la o anumita scara reprezinta curentul de colector, este distorsionata Acest lucru se produce deoarece pentru circuitul neliniar cu tranzistor se poate folosi modelul liniar in curent alternativ numai pentru functionarea la semnale mici E2=20 0 V corespunde functionarii la semnale mari 3 INSTRUCȚIUNI SPICE Se folosesc urmatoarele linii: a) Semnalul sinusoidal amortizat este descris astfel: Pentru sursele independente forma generala este: V SIN (voff vampl freq td df phase) I SIN (ioff iampl freq td df phase) I= ioff + iampl*expi-(t-td)*df*sin(2*3 14*freq*(t-td)+3 14*Phase/180) unde td [ [ ]][SKIPBP] =pasul de timp pentru tiparirea rezultatelor; = timpul final; = timpul de la care se incepe tiparirea (vizualizarea) rezultatelor; 77 =limita superioara a pasului de timp; daca nu se specifica se considera ceiling=ftime/50 pentru un circuit care conține elemente dinamice sau ceiling=pstep pentru un circuit rezistiv SKIPBP = utilizarea conditiilor initiale impuse (U= pentru condensator si I= pentru bobina, care sunt specificate in liniile care definesc elementele respective) c) definirea unui subcircuit SUBCKT [ *] [PARAMS: [= ]* ] (descrierea elementelor) ENDS Utilizarea subcircuitului : X [ ]* [PARAMS: = *>] = numele modelului (definit in linia SUBCKT) = numele elementului modelat cu subcircuitul respectiv; = nodurile prin care subcircuitul se conecteaza cu circuitul 78 LUCRAREA 6 REGIMUL TRANZITORIU AL CIRCUITELOR NELINIARE 1 PROBLEME 1 ) Pentru circuitul din figura sa se faca analiza Fourier a tensiunii de iesire pentru E2 =1 mV si E2 =10 V incazul in care frecventa sursei e2 este f=1 kHz Rb = 20 kQ ;C=100 uF, R, = 100 Q; E1 = 5V; Vcc = 9V; e2 = E2 sin (rot); TB: в= 200 ; Cje = 200 pF; Cjc = 200 pF; Vje = 0 6V; Vjc = 0 6V; 2 ) Sa se simuleze convertorul dc - dc : L=20mH ; C=1 25uF; R=400Q ; E=100V; Comutator: - frecventa de comutatie f= 10kHz ; - durata de inchidere reprezinta 50% din perioada de comutatie; - rezistente : Rinchis = 1Q ; Rdeschis = 1e+06Q ; D - dioda avand modelul predefinit de PSPICE; Sa se vizualizeze tensiunea de comanda a comutatorului(pe un interval de 5 perioade de comutatie) si tensiunea la iesirea convertorului(atat pe intreaga durata a regimului tranzitoriu cat si 3 perioade din regimul permanent) Sa se traseze spectrul de frecvente pentru tensiunea de comanda a comutatorului si pentru tensiunea la iesirea convertorului 2 INTERPRETAREA REZULTATELOR 79 La problema 1 , realizând analiza Fourier a tensiunii de iesire rezulta fig 1 pentru E2=1 0 mV si fig 2 pentru E2=10 0 V □ V(5)- V(4) Frequency o V(5,4) Frequency Fig 2 La problema 2 chestiunile care se vor urmari sunt: - modul de utilizare a comutatorului; - descrierea semnalelor de tip puls pentru surse in regim tranzitoriu; -trasarea spectrului de frecvente al unei marimi periodice 80 110 0 mV o V(4) Ti me Fig 3 In fig 3 se prezintă semnalul generat cu funcția pulse, iar in fig 4 se prezintă semnalul la ieșirea convertorului Spectrele tensiunii de comanda si pentru tensiunea de iesire sunt prezentate in fig 5 - 7 Figura 7 contine un detaliu al spectrului tensiunii de iesire care permite evidentierea armonicelor frecventei de comanda x V(3) Ti me Fig 4 81 8 0 0mV 4 0 0mV 0V 0Hz o V(4) 10 0KH z 2 0 0 KH z Freque ncy 3 0 0 KH z 4 0 0 KHz Fig 5 Se observa ca spectrul tensiunii de comanda este discret (conține numai componentele armonice ale frecventei de comutatie si componenta de current continuu) deoarece aceasta marime este periodica Tensiunea la iesirea convertorului nefiind periodica (deoarece contine si componente tranzitorii care nu s-au amortizat) spectrul acesteia nu este discret o V(3) Frequency Fig 6 Pentru trasarea spectrului tensiunii de iesire s-au considerat numai ultimile 5 perioade din raspuns (o forma de unda aproapiata de o functie periodica) 82 Frequeno у Fig 7 3 INSTRUCȚIUNI SPICE Se vor utiliza următoarele linii: a) Semnalul de tip puls pentru surse in regim tranzitoriu este descris de: PULSEfii i2 td trise tfaii pw per) pentru sursa de curent PULSE(vi v2 td trise tfaii pw per) pentru sursa de tensiune td: per : C C > Mărime Semnificație Valoare predefinita Unitate de masura h Valoare inițiala - A іг Valoare de vârf - A ta Timp de intarziere 0 0 s trise Timp de creștere TSTEP s tfall Timp de descreștere TSTEP s pw Durata impulsului TSTOP s per Perioada TSTOP s 83 b) comutatorul comandat in tensiune S = nodurile intre care se conectează; = nodurile tensiunii de comanda; = numele modelului; MODEL VSWITCH ( RON= , ROFF= , VON= , VOFF= ) Valori predefinite : RON IQ ; ROIT 1V1 * = frecventa considerata fundamentala ; = semnalele pentru care se face analiza Fourier; 84 LUCRAREA 7 COMPORTAREA CALITATIVA A CIRCUITELOR DE ORDINUL II 1 PROBLEME 1 ) Sa se simuleze circuitul cu condițiile inițiale si valorile din tabelul 1 : L=1H; C=1F; Caz R [Q] Condiții inițiale (ис(0),і[(0)) 1 3 (1,1) (-1,1) (1,-1) (-1,-1) 2 -3 (1,1) (-1,1) (1,-1) (-1,-1) 3 1 (1,1) (-1,1) (1,-1) (-1,-1) 4 -1 (1,1) (-1,1) (1,-1) (-1,-1) 5 le-9 (1,1) (-1,1) (1,-1) (-1,-1) Tabelul 1 Sa se indice in cazurile 1 5 tipul comportării calitative a circuitului in jurul punctului de echilibru 2 a) Sa se simuleze circuitul din figura pentru i (0) = 6A Sa se vizualizeze i(t) si curba i = f(V( 1)) si sa se explice rezultatul obtinut 0 b) Intre nodurile 1 si 0 se conectează un condensator cu C= ImF avand condiția inițiala Uc0= 8V Sa se vizualizeze traiectoria in planul fazelor (iL - uc) Sa se compare rezultatele obtinute pentru eroarea relative RELTOL cu valoarea predefinita( le-05) si valoarea impusa de le-07 2 INTERPRETAREA REZULTATELOR 85 1 Circuitul din problema 1 este un circuit liniar, dinamic, de ordinul doi, cu excitație nula, cu ecuațiile: Uq + Ri + UL = 0 ?c = iL = i Tinand cont de ecuatiile de functionare ale bobinei si condensatorului UL = LiL ' ic = CUc rezulta ecuatiile de stare ale circuitului: "Uq iL Uq iL Pentru a studia comportarea calitativa a unui circuit liniar cu excitatia nula se determina valorile proprii s1 si s2 ale matricei A a ecuatiei X = Ax Valorile proprii s1 si s2 sunt solutiile ecuatiei: det C R 5 L respectiv care pentru L=1H si C=1F devine s2 + Rs +1 = 0 Punctele de echilibru se determina rezolvand sistemul de ecuatii Ax=0 adica a11 a21 X1Q a12 a22 J[ x2Q In cazul problemei A = -1 si A = 1 deci A ф 0 si sistemul admite numai solutia banala si 0 - 5 L 0 1 - R 0 1 L C R L 1 = 0 2 R s + — L 5 + —^ = 0 L ■ C = 0 originea este singurul punct de echilibru adica x1q = x2q = 0 Evolutia circuitului plecand de la o stare initiala(uC(0) si iL(0)) poate fi reprezentata printr-o curba numita traiectorie de faza in planul de coordonate uC si iL numit planul fazelor Evolutia circuitului corespunzatoare mai multor stari initiale poate fi reprezentataprintr-o multime de traiectorii in planul fazelor Aceste traiectorii formeazaun portret de faza Considerand perechile de valori UC0 si Il0 din tabelul 1 se determinatraiectoriile corespunzatoare starilor initiale date, care formeaza portrete de faza in cele 4 cazuri considerate In functie de valorile rezistentei R din tabelul 1 rezulta 4 cazuri Cazul1 R = 3 Q si valorile proprii sunt: -3+45 -3-45 s1 = - s2 = - 1 2 2 2 deci s2 s1 > 0 si originea este un nod instabil Cazul 3 R = 1Q si valorile proprii sunt: 1 45 1 45 S1 = -+ j - S2 = j 1 2 2 2 2 2 deci constanta de atenuare a = - у si pentru a 0, punctul de echilibru este un focar instabil 87 Soluția x(t) a ecuatiei de stare X = Ax, are doua componente: uC(t) si iL(t) Cazul 5 R~0 si valorile proprii sunt 2 = ± j si originea este centru 2 Pentru circuitul dinamic neliniar de la problema 2 punctul a), parcursul dinamic este: unde punctele Q1(1,-1) si Q2(-1,1) sunt puncte de impas iar originea 0 este un nod instabil Simuland circuitul cu PSPICE se constata ca traiectoria evolueaza pe curba plecand din punctul corespunzator conditiei initiale iL(0) si se opreste in punctul de impas intalnit Q1(1,-1) La punctul b) al problemei 2 se constata ca legând la bornele circuitului un condensator C, circuitul se comporta ca un oscilator de relaxare In cazul circuitului real nu este necesara introducerea in circuit a condensatorului C, rolul acestuia fiind indeplinit de capacitatea parazita dintre bornele circuitului Cel de-al doilea punct al problemei 2 are ca scop ilustrarea preciziei analizei efectuata cu PSPICE in functie de diferite valori date pentru eroarea relativa RELTOL 88 RELTOL=1e-05 RELTOL=1e-07 3 INSTRUCȚIUNI SPICE OPTIONS [ *] [ = *] - setarea unor opțiuni: Flag Options ACCT summary & accounting EXPAND show subcircuit expansion LIBRARY list lines from library files LIST output summary NODE output netlist NOECHO suppress listing NOMOD suppress model parameter listing NOPAGE suppress banners OPTS output option values Value Options ABSTOL best accuracy of currents CHGTOL best accuracy of charges CPTIME CPU time allowed DEFAD MOSFET default AD DEFAS MOSFET default AS DEFL MOSFET default L 89 DEFW MOSFET default W GMIN minimal conductance, any branch ITL1 iteration number - DC & bias point blind limit ITL2 iteration number - DC & bias point guess limit ITL4 iteration number - transient per-point limit ITL5 iteration number - transient total, all points LIMPTS maximal number of points for print/plot NUMDGT number of digits in output file PIVREL relative magnitude for matrix pivot PIVTOL absolute magnitude for matrix pivot RELTOL relative accuracy of voltages and currents TNOM default temperature TRTOL transient accuracy adjustment VNTOL best accuracy of voltages WIDTH printing width in output file 90 LUCRAREA 8 PROPRIETATI CALITATIVE IN CIRCUITELE ELECTRICE NELINIARE IN REGIM TRANZITORIU 1 PROBLEME Sa se simuleze circuitul cu următorii parametri: - dioda D : IS=8 3e-15A; RS=9 6Q; TT=4e-06s (timp de tranzit) ; CJO=300pF ; M=0 4 (gradient); VJ=0 75V (potential de jonctiune); R = 15Q ; L = 10mH; e(t)= E0sin (cot) ; f=100kHz o 3 Se vor utiliza un pas de afisare de T/2000=5 e-09 s si o valoare RELTOL=1e-6 Sa se reprezinte grafic dependenta i(t) = f ( e(t)) si sa se interpreteze rezultatele in urmatoarele situatii: a) E0 = 0 5 V Se vor simula 40 perioade , dintre care se vor afisa ultimele 4 b) E0 = 1 5 V Se vor simula 40 perioade, dintre care se vor afisa ultimele 4 c) E0 = 6 5 V Se vor simula 60 de perioade, dintre care se vor afisa ultimele 20 2 INTERPRETAREA REZULTATELOR In urma simularii se determina raspunsul I(L1) pentru cele trei cazuri: a) pentru E0=0,5 V rezulta un raspuns periodic de frecventa excitatiei 91 4 0 0uA- 400uA' -4 0 0uA-| 1 -1 -1 - 360us 3 70us 380us 390us 400us □ I(l1) T i m e -400uA - -500mV 0V 500mV o I(l1) V ( 1 ) Fig 1 Eo=O,5 V Comportarea circuitului in acest caz este o comportare obișnuita deoarece pentru amplitudini relativ mici ale excitatiei comportarea diodei poate fi aproximata prin modelul de semnal mic b) pentru Eo=1,5 V se observa raspunsul subarmonic(periodic de frecventa W2) O I(l1) Ti me Fig 2 Eo=1,5 V V( 1 ) c) pentru E0=6,5 V se observa comportarea haotica a circuitului(raspuns neperiodic la excitatie periodica) 92 Fig 3 Eo=6,5 V □ I(l1) T i m e o I (11) V( 1 ) 93 CAP 4 SOLUȚIILE PROBLEMELOR LUCRAREA 1 Problema 1 * * * * CIRCUIT DESCRIPTION *circuit de CC Ll/pr l R1 1 23 R24 0 6 R4 3 2 6 R5 4 5 3 13 42 1 E4 3 5 4 2 9 VI 1 09 V6 5 0 12 OP DC VI 0 10 1 PROBE END **** SMALL SIGNAL BIAS SOLUTION TEMPERATURE = 27 000 DEG C NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE ( 1) 9 0000 ( 2) 7 5000 ( 3) -1 5000 ( 4) 6 0000 ( 5) 12 0000 VOLTAGE SOURCE CURRENTS NAME CURRENT VI -5 000E-01 V6 -5 000E-01 TOTAL POWER DISSIPATION 1 05E+01 WATTS **** VOLTAGE-CONTROLLED VOLTAGE SOURCES NAME E4 V-SOURCE -1 350E+01 I-SOURCE 1 500E+00 JOB CONCLUDED 94 Problema 2 pct a * * * * CIRCUIT DESCRIPTION *circuit de CC Ll/pr 2 R1 0 1 1 R2 2 3 1 R3 3 02 II 02 1 VI 2 1 2 HI 4 3 VI -2 V44 0 1 *R4 4 0 2 OP END **** SMALL SIGNAL BIAS SOLUTION TEMPERATURE = 27 000 DEG C NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE ( 1) -2 0000 ( 2) 0 0000 ( 3) -3 0000 ( 4) 1 0000 VOLTAGE SOURCE CURRENTS NAME CURRENT VI -2 000E+00 V4 4 500E+00 TOTAL POWER DISSIPATION -5 00E-01 WATTS **** OPERATING POINT INFORMATION TEMPERATURE = 27 000 DEG C **** CURRENT-CONTROLLED VOLTAGE SOURCES NAME HI V-SOURCE 4 000E+00 I-SOURCE -4 500E+00 pct b * * * * CIRCUIT DESCRIPTION *circuit de CC Ll/pr 2 R1 0 1 1 R2 2 3 1 R3 3 02 II 02 1 VI 2 1 2 95 HI 4 3 VI -2 V44 0 2 *R4 4 О 2 OP END ERROR — Convergence problem in bias point calculation Lasl node voltages tried were: NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE ( 1) 10 00E+09 ( 2) 10 00E+09 ( 3) 10 00E+09 ( 4) 2 0000 These voltages failed to converge: V(l) = le+010 / le+010 V(2) = le+010/le+010 V(3) = le+010/le+010 These supply currents failed to converge: I(V1) = le+010/le+010 I(V4) = -le+010/-le+010 I(H1) = le+010/le+010 Try running with OPT1ON STEPGMIN JOB CONCLUDED TOTAL JOB TIME 01 pct c * * * * CIRCUIT DESCR1PT1ON *circuit de ce Ll/pr 2 R1 0 1 1 R22 3 1 R3 3 0 2 11 02 1 VI 2 1 2 HI 4 3 VI -2 *V4 4 0 2 R4402 OP END **** SMALL SIGNAL BIAS SOLUT1ON TEMPERATURE = 27 000 DEG C 96 NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE ( l)-2 000E-12 ( 2) 2 0000 ( 3) 1 0000 ( 4) 1 0000 VOLTAGE SOURCE CURRENTS NAME CURRENT VI -2 000E-12 TOTAL POWER DISSIPATION 4 00E-12 WATTS **** CURRENT-CONTROLLED VOLTAGE SOURCES NAME HI V-SOURCE 4 000E-12 I-SOURCE -5 000E-01 JOB CONCLUDED LUCRARAEA 2 Problema 1 * * * * CIRCUIT DESCRIPTION *Circ cu dioda Zener cdzl cir R1 1 2 Ik R2 2 0 2k D 0 2 Dioda Zener *Modelul diodei *Bv=4 tens de strapungere inversa MODEL Dioda Zener D(Is=10f Rs=0 + Bv=4 Ibv= 10m Cj o=0) *Tensiunea de polarizare V1 1 0 1 *Analiza in curent continuu DC LIN V1 -10 15 1 PROBE END **** Diode MODEL PARAMETERS Dioda Zener IS 10 000000E-15 BV 4 IBV 01 97 V( 1 ) Problema 2 pct a * * * * CIRCUIT DESCRIPTION *circuit cu rezistenta neliniara cnll cir R1 1 2 3 V1 1 0 10 G1 2 0 TABLE {V(2)}= +(0,0), (2,5), (4,0 5), (8,9) NODESET V(2)=0 0 OP END **** SMALL SIGNAL BIAS SOLUTION TEMPERATURE = 27 000 DEG C NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE ( 1) 10 0000 ( 2) 1 1765 VOLTAGE SOURCE CURRENTS NAME CURRENT V1 -2 941E+00 TOTAL POWER DISSIPATION 2 94E+01 WATTS **** VOLTAGE-CONTROLLED CURRENT SOURCES NAME G1 I-SOURCE 2 941E+00 98 pct b * * * * CIRCUIT DESCRIPTION *circuit cu rezistenta neliniara cnll cir R1 1 2 3 V1 1 0 10 G1 2 0 TABLE {V(2)}= +(0,0), (2,5), (4,0 5), (8,9) NODESET V(2)=3 0 OP END **** SMALL SIGNAL BIAS SOLUTION TEMPERATURE = 27 000 DEG C NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE ( 1) 10 0000 ( 2) 3 2174 VOLTAGE SOURCE CURRENTS NAME CURRENT V1 -2 261E+00 TOTAL POWER DISSIPATION 2 26E+01 WATTS **** VOLTAGE-CONTROLLED CURRENT SOURCES NAME G1 I-SOURCE 2 261E+00 pct c * * * * CIRCUIT DESCRIPTION *circuit cu rezistenta neliniara cnl1 cir R1 1 2 3 V1 1 0 10 G1 2 0 TABLE {V(2)}= +(0,0), (2,5), (4,0 5), (8,9) NODESET V(2)=7 0 OP END **** SMALL SIGNAL BIAS SOLUTION TEMPERATURE = 27 000 DEG C 99 NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE ( 1) 10 0000 ( 2) 4 6102 VOLTAGE SOURCE CURRENTS NAME CURRENT V1 -1 797E+00 TOTAL POWER DISSIPATION 1 80E+01 WATTS **** VOLTAGE-CONTROLLED CURRENT SOURCES NAME G1 I-SOURCE 1 797E+00 pct d * * * * CIRCUIT DESCRIPTION *circuit cu rezistenta neliniara cnl1 cir R1 1 2 3 V1 1 0 10 G1 2 0 TABLE {V(2)}= +(0,0), (2,5), (4,0 5), (8,9) * NODESET V(2)=7 0 OP END **** SMALL SIGNAL BIAS SOLUTION TEMPERATURE = 27 000 DEG C NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE ( 1) 10 0000 ( 2) 1 1765 VOLTAGE SOURCE CURRENTS NAME CURRENT V1 -2 941E+00 TOTAL POWER DISSIPATION 2 94E+01 WATTS **** VOLTAGE-CONTROLLED CURRENT SOURCES NAME G1 I-SOURCE 2 941E+00 100 LUCRAREA 3 Punctul a **** CIRCUIT DESCRIPTION *etb31 cir circ cu tranz bipolare * options reltol=1e-6 itl5=0 numdgt=8 V1 2 0 DC 5 V2 5 0 DC 9 R1 2 3 1MEG R2 5 4 4K Q1 4 3 0 TB MODEL TB NPN(Is=10 f Bf=200 Cje=2p + Cj c=2p Vje= 6 Vj c= 6) DC LIN V1 0 25 05 * PRINT DC V(R2) PROBE END **** BJT MODEL PARAMETERS TB NPN IS 10 000000E-15 BF 200 NF 1 BR 1 NR 1 CJE 2 000000E-12 VJE 6 CJC 2 000000E-12 VJC 6 CN 2 42 D 87 JOB CONCLUDED 101 V1 Punctul b **** CIRCUIT DESCRIPTION *etb32 cir circ cu tranz bipolare V1 2 1 DC 5 V 1 0 + DC 0 + AC 1 0 V2 5 0 DC 9 R1 2 3 1MEG R2 5 4 4K Q1 4 3 0 TB MODEL TB NPN(Is=10 f Bf=200 Cje=2p + Cjc=2p Vje= 6 Vjc= 6) * DC LIN V1 0 25 05 AC DEC 100 1 100MEG OP PROBE END **** BJT MODEL PARAMETERS TB NPN IS 10 000000E-15 BF 200 NF 1 BR 1 NR 1 CJE 2 000000E-12 VJE 6 CJC 2 000000E-12 VJC 6 CN 2 42 102 D 87 **** SMALL SIGNAL BIAS SOLUTION TEMPERATURE NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE ( 1) 0 0000 ( 2) 5 0000 ( 3) 6515 ( 4) 5 5212 ( 5) 9 0000 27 000 DEG C NODE VOLTAGE VOLTAGE SOURCE CURRENTS NAME CURRENT V1 -4 349E-06 V -4 349E-06 V2 -8 697E-04 TOTAL POWER DISSIPATION 7 85E-03 WATTS **** BIPOLAR JUNCTION TRANSISTORS NAME Q1 MODEL TB IB 4 35E-06 IC 8 70E-04 VBE 6 51E-01 VBC -4 87E+00 VCE 5 52E+00 BETADC 2 00E+02 GM 3 36E-02 RPI 5 95E+03 RX 0 00E+00 RO 1 00E+12 CBE 3 49E-12 CBC 9 64E-13 CJS 0 00E+00 BETAAC 2 00E+02 CBX/CBX2 0 00E+00 FT/FT2 1 20E+09 103 8 0 0mV 4 0 0mV 0V 100mHz 1 0Hz □ V(5)- V(4) 100Hz 10KHz 1 0MHz Frequency 100MHz **** CIRCUIT DESCRIPTION *etb32 cir circ cu tranz bipolare V1 2 1 DC 5 V 1 0 + DC 0 + AC 20 0 V2 5 0 DC 9 R1 2 3 1MEG R2 5 4 4K Q1 4 3 0 TB MODEL TB NPN(Is=10 f Bf=200 Cje=2p + Cj c=2p Vje= 6 Vj c= 6) * DC LIN V1 0 25 05 AC DEC 100 1 100MEG OP PROBE END 1 1 1 2 000000E-12 6 2 000000E-12 6 2 42 87 **** BJT MODEL PARAMETERS TB NPN IS 10 000000E-15 BF 200 NF BR NR CJE VJE CJC VJC CN D 104 **** SMALL SIGNAL BIAS SOLUTION TEMPERATURE NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE ( 1) 0 0000 ( 2) 5 0000 ( 3) 6515 ( 4) 5 5212 ( 5) 9 0000 27 000 DEG C NODE VOLTAGE VOLTAGE SOURCE CURRENTS NAME CURRENT V1 -4 349E-06 V -4 349E-06 V2 -8 697E-04 TOTAL POWER DISSIPATION 7 85E-03 WATTS **** BIPOLAR JUNCTION TRANSISTORS NAME Q1 MODEL TB IB 4 35E-06 IC 8 70E-04 VBE 6 51E-01 VBC -4 87E+00 VCE 5 52E+00 BETADC 2 00E+02 GM 3 36E-02 RPI 5 95E+03 RX 0 00E+00 RO 1 00E+12 CBE 3 49E-12 CBC 9 64E-13 CJS 0 00E+00 BETAAC 2 00E+02 CBX/CBX2 0 00E+00 FT/FT2 1 20E+09 105 20V 10V 0V 100mHz 1 0Hz □ V(5)- V(4) 100Hz 10KHz 1 0MHz 100MHz Frequency LUCRARAEA 4 Problema 1 **** CIRCUIT DESCRIPTION *circ lin c a CLA1 CIR R1 2 3 10 R2 3 5 1 R3 3 4 2 L1 1 2 50m L2 6 7 25m L3 0 4 50m K1 l1 l3 0 5 K2 l2 l3 0 707 C2 5 6 1 5m V SIN 7 0 +DC 0 +AC 30 90 I SIN 0 1 +DC 0 +AC 2 0 AC LIN 1 50 50 PRINT AC IM(R2) IP(R2) IR(R2) II(R2) PROBE END NODE VOLTAGE NODE **** SMALL SIGNAL BIAS SOLUTION TEMPERATURE = 27 000 DEG C VOLTAGE NODE VOLTAGE NODE VOLTAGE ( 1) 0 0000 ( 2) ( 5) 0 0000 ( 6) 0 0000 0 0000 ( 3) ( 7) 0 0000 ( 4) 0 0000 0 0000 VOLTAGE SOURCE CURRENTS 106 NAME CURRENT V SIN 0 000E+00 TOTAL POWER DISSIPATION 0 00E+00 WATTS **** AC ANALYSIS TEMPERATURE = 27 000 DEG C FREQ IM(R2) IP(R2) IR(R2) II(R2) 5 000E+01 1 138E-01 -6 592E+01 4 645E-02 -1 039E-01 Problema 2a **** CIRCUIT DESCRIPTION *CIRC LIN DE CA FIL1 R1 1 2 1K L1 2 3 3M L2 3 4 12M C1 3 0 8N C2 4 0 11 5N V1 1 0 +DC 0 +AC 1 0 +SIN(0 1 50 0 0 0) AC LIN 100 1 1MEG PROBE END **** SMALL SIGNAL BIAS SOLUTION TEMPERATURE = 27 000 DEG C NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE ( 1) 0 0000 ( 2) 0 0000 ( 3) 0 0000 ( 4) 0 0000 VOLTAGE SOURCE CURRENTS NAME CURRENT V1 0 000E+00 TOTAL POWER DISSIPATION 0 00E+00 WATTS 107 Frequency Problema 2b **** CIRCUIT DESCRIPTION *CIRC LIN DE CA FIL2 R1 1 2 1 5K L1 3 0 2M L2 4 0 12M C1 2 3 10N C2 3 4 8U V1 1 0 +DC 0 +AC 1 0 +SIN(0 1 50 0 0 0) AC LIN 100 1 2MEG PROBE END **** SMALL SIGNAL BIAS SOLUTION TEMPERATURE = 27 000 DEG C NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE ( 1) 0 0000 ( 2) 0 0000 ( 3) 0 0000 ( 4) 0 0000 VOLTAGE SOURCE CURRENTS NAME CURRENT V1 0 000E+00 TOTAL POWER DISSIPATION 0 00E+00 WATTS 108 1 0V- 0 5V- 0V- 1 0Hz □ V(4) 100Hz 10KHz 1 0MHz 100MHZ Frequency Problema 2c **** CIRCUIT DESCRIPTION *CIRC LIN DE CA FIL3 R1 1 2 1 5K RS 4 0 1 5K L1 2 3 2M C1 2 0 10N C2 3 4 8U C3 4 0 1N V1 1 0 +DC 0 +AC 1 0 +SIN(0 1 50 0 0 0) AC DEC 100 1 2MEG PROBE END **** SMALL SIGNAL BIAS SOLUTION TEMPERATURE = 27 000 DEG C NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE ( 1) 0 0000 ( 2) 0 0000 ( 3) 0 0000 ( 4) 0 0000 VOLTAGE SOURCE CURRENTS NAME CURRENT V1 0 000E+00 TOTAL POWER DISSIPATION 0 00E+00 WATTS 109 □ V(4) Frequency Problema 2d **** CIRCUIT DESCRIPTION *CIRC LIN DE CA FIL4 R1 1 2 1 5K L1 5 0 10M L2 2 3 20M L3 4 0 4M C1 2 5 1U C2 2 3 40U C3 3 4 1U V1 1 0 +DC 0 +AC 1 0 +SIN(0 1 50 0 0 0) AC DEC 100 1 2MEG PROBE END **** SMALL SIGNAL BIAS SOLUTION TEMPERATURE = 27 000 DEG C NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE ( 1) 0 0000 ( 2) 0 0000 ( 3) 0 0000 ( 4) 0 0000 ( 5) 0 0000 VOLTAGE SOURCE CURRENTS NAME CURRENT V1 0 000E+00 110 TOTAL POWER DISSIPATION 0 00E+00 WATTS O V(3) Frequency LUCRAREA 5 Problema 1 **** CIRCUIT DESCRIPTION *est51 cir CIR cu AO/ SUBcirc V1 1 0 + DC 0 + AC 1m 0 + SIN(0 1m 500 0 0 0) * R1 1 2 1 C1 1 2 1 R2 2 3 1 C2 3 6 1 R3 2 4 1 R4 4 5 1 R5 5 6 1 C5 5 6 1 R6 1 7 1 C6 7 5 1 * X1 AO X2 AO 111 * SUBCKT AO 2 1 4 Ri 1 2 1MEG Re 3 4 1m Ce 4 2 1p E 3 2 TABLE { V(2)-V(1)} = (- 0015,-15)(0,0)( 0015,15) ENDS * AC DEC 100 1 10MEG * tran 1n 4 m probe options reltol=1e-6 itl5=0 numdgt=8 end **** SMALL SIGNAL BIAS SOLUTION TEMPERATURE = 27 000 DEG C NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE ( 1) 0 00000000 ( 2) 0 00000000 ( 3) 0 00000000 ( 4) 0 00000000 ( 5) 0 00000000 ( 6) 0 00000000 ( 7) 0 00000000 ( X1 3) 0 00000000 ( X2 3) 0 00000000 VOLTAGE SOURCE CURRENTS NAME CURRENT V1 0 000E+00 TOTAL POWER DISSIPATION 0 00E+00 WATTS □ V(6) Frequency 112 1 0 0 5 0-| -— 100mHz 1 0Hz 100Hz 10KHz 1 0MHz o V(6)/ V(1) Frequency *est51 cir CIR cu AO/ SUBcirc V1 1 0 + DC 0 + AC 1m 0 + SIN(0 1m 500 0 0 0) * R1 1 2 1 C1 1 2 1 R2 2 3 1 C2 3 6 1 R3 2 4 1 R4 4 5 1 R5 5 6 1 C5 5 6 1 R6 1 7 1 C6 7 5 1 * X1 AO X2 AO * SUBCKT AO 2 1 4 Ri 1 2 1MEG Re 3 4 1m Ce 4 2 1p E 3 2 TABLE { V(2)-V(1)} = (- 0015,-15)(0,0)( 0015,15) ENDS * * AC DEC 100 1 10MEG tran 1n 4 m probe options reltol=1e-6 itl5=0 numdgt=8 end **** INITIAL TRANSIENT SOLUTION TEMPERATURE = 27 000 DEG C 113 NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE ( 1) 0 00000000 ( 2) 0 00000000 ( 3) 0 00000000 ( 4) 0 00000000 ( 5) 0 00000000 ( б) 0 00000000 ( 7) 0 00000000 ( X1 3) 0 00000000 ( X2 3) 0 00000000 VOLTAGE SOURCE CURRENTS NAME CURRENT V1 0 000E+00 TOTAL POWER DISSIPATION 0 00E+00 WATTS o V(6) Ti me **** CIRCUIT DESCRIPTION *est51 cir CIR cu AO/ SUBcirc V1 1 0 + DC 0 + AC 1m 0 + SIN(0 10 500 0 0 0) * R1 1 2 1 C1 1 2 1 R2 2 3 1 C2 3 б 1 R3 2 4 1 R4 4 5 1 R5 5 б 1 C5 5 б 1 R6 1 7 1 C6 7 5 1 * X1 AO 114 X2 AO * SUBCKT AO 2 1 4 Ri 1 2 1MEG Re 3 4 1m Ce 4 2 1p E 3 2 TABLE { V(2)-V(1)} = (- 0015,-15)(0,0)( 0015,15) ENDS * * AC DEC 100 1 10MEG tran 1n 4 m probe options reltol=1e-6 itl5=0 numdgt=8 end **** INITIAL TRANSIENT SOLUTION TEMPERATURE = 27 000 DEG C NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE ( 1) 0 00000000 ( 2) 0 00000000 ( 3) 0 00000000 ( 4) 0 00000000 ( 5) 0 00000000 ( 6) 0 00000000 ( 7) 0 00000000 ( X1 3) 0 00000000 ( X2 3) 0 00000000 VOLTAGE SOURCE CURRENTS NAME CURRENT V1 0 000E+00 o V(6) Ti me Problema 2 **** CIRCUIT DESCRIPTION 115 *etb52 cir circ cu tranz bipolare options reltol=1e-6 V1 1 0 DC 5 C3 2 3 100u *V3 3 0 AC 1m SIN(0 1m 1K 0 0) V3 3 0 AC 10 SIN(0 10 1K 0 0) V2 5 0 DC 9 R1 1 2 20K R2 5 4 100 Q1 4 2 0 TB MODEL TB NPN(IS=10 F BF=200 CJE=200p + CJC=200p) * DC LIN V1 0 25 05 * AC DEC 100 1 100MEG TRAN 1u 10m OP PROBE END **** BJT MODEL PARAMETERS TB NPN IS 10 000000E-15 BF 200 NF 1 BR 1 NR 1 CJE 200 000000E-12 CJC 200 000000E-12 CN 2 42 D 87 **** SMALL SIGNAL BIAS SOLUTION TEMPERATURE = 27 000 DEG C NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE ( 1) 5 0000 ( 2) 7521 ( 3) 0 0000 ( 4) 4 7521 ( 5) 9 0000 VOLTAGE SOURCE CURRENTS NAME CURRENT V1 -2 124E-04 116 V3 0 000E+00 V2 -4 248E-02 TOTAL POWER DISSIPATION 3 83E-01 WATTS **** OPERATING POINT INFORMATION TEMPERATURE = 27 000 DEG C **** BIPOLAR JUNCTION TRANSISTORS NAME Q1 MODEL TB IB 2 12E-04 IC 4 25E-02 VBE 7 52E-01 VBC -4 00E+00 VCE 4 75E+00 BETADC 2 00E+02 GM 1 64E+00 RPI 1 22E+02 RX 0 00E+00 RO 1 00E+12 CBE 3 35E-10 CBC 1 09E-10 CJS 0 00E+00 BETAAC 2 00E+02 CBX/CBX2 0 00E+00 FT/FT2 5 89E+08 **** INITIAL TRANSIENT SOLUTION TEMPERATURE = 27 000 DEG C NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE ( 1) 5 0000 ( 2) 7521 ( 3) 0 0000 ( 4) 4 7521 ( 5) 9 0000 VOLTAGE SOURCE CURRENTS NAME CURRENT V1 -2 124E-04 V3 0 000E+00 V2 -4 248E-02 TOTAL POWER DISSIPATION 3 83E-01 WATTS 117 4 50V 4 00V 4 25V 0 s 2ms 4ms 6ms 8ms 10ms □ V(5)- V(4) Ti me □ V(5)- V(4) Ti me LUCRAREA 6 Problema 1 **** CIRCUIT DESCRIPTION *etb61 cir circ cu tranz bipolare options reltol=1e-6 V1 1 0 DC 5 118 C3 2 3 100u V3 3 0 AC 1m SIN(0 1m 1K 0 0) *V3 3 0 AC 10 SIN(0 10 1K 0 0) V2 5 0 DC 9 R1 1 2 20K R2 5 4 100 Q1 4 2 0 TB MODEL TB NPN(IS=10 F BF=200 CJE=200p + CJC=200p) * TRAN 1u 10m TRAN 1u 10m 5m FOUR 1K 10 V(5,4) * OP PROBE END **** BJT MODEL PARAMETERS TB NPN IS 10 000000E-15 BF 200 NF 1 BR 1 NR 1 CJE 200 000000E-12 CJC 200 000000E-12 CN 2 42 D 87 **** INITIAL TRANSIENT SOLUTION TEMPERATURE = 27 000 DEG C NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE ( 1) 5 0000 ( 2) 7521 ( 3) 0 0000 ( 4) 4 7521 ( 5) 9 0000 VOLTAGE SOURCE CURRENTS NAME CURRENT V1 -2 124E-04 V3 0 000E+00 V2 -4 248E-02 TOTAL POWER DISSIPATION 3 83E-01 WATTS **** FOURIER ANALYSIS TEMPERATURE = 27 000 DEG C 119 FOURIER COMPONENTS OF TRANSIENT RESPONSE V(5,4) DC COMPONENT = 4 247607E+00 HARMONIC FREQUENCY FOURIER NORMALIZED PHASE NORMALIZED NO (HZ) COMPONENT COMPONENT (DEG) PHASE (DEG) 1 1 000E+03 1 625E-01 1 000E+00 7 960E-01 0 000E+00 2 2 000E+03 1 589E-03 9 777E-03 -8 396E+01 -8 476E+01 3 3 000E+03 5 077E-04 3 125E-03 -5 091E+01 -5 171E+01 4 4 000E+03 1 470E-04 9 045E-04 1 237E+02 1 229E+02 5 5 000E+03 3 713E-04 2 285E-03 7 320E+01 7 240E+01 6 6 000E+03 1 359E-04 8 364E-04 -1 204E+02 -1 212E+02 7 7 000E+03 2 854E-04 1 756E-03 -1 695E+02 -1 703E+02 8 8 000E+03 1 450E-04 8 922E-04 -4 044E+00 -4 840E+00 9 9 000E+03 2 411E-04 1 484E-03 -4 427E+01 -4 506E+01 10 1 000E+04 1 767E-04 1 088E-03 1 118E+02 1 110E+02 TOTAL HARMONIC DISTORTION = 1 092516E+00 PERCENT 5 0V 2 5V 0V ] A 0Hz 5KHz □ V(5)- V(4) 10KHz 15KHz Frequency 20KHz 25KHz 30KHz 120 1 0V 0 5V 0Hz 2KHz □ V(5)- V(4) 4KHz 6KHz Frequency 8KHz 10KHz **** CIRCUIT DESCRIPTION *etb61 cir circ cu tranz bipolare options reltol=1e-6 V1 1 0 DC 5 C3 2 3 100u *V3 3 0 AC 1m SIN(0 1m 1K 0 0) V3 3 0 AC 10 SIN(0 10 1K 0 0) V2 5 0 DC 9 R1 1 2 20K R2 5 4 100 Q1 4 2 0 TB MODEL TB NPN(IS=10 F BF=200 CJE=200p + CJC=200p) * TRAN 1u 10m TRAN 1u 10m 5m FOUR 1K 10 V(5,4) * OP PROBE END **** BJT MODEL PARAMETERS TB NPN IS 10 000000E-15 BF 200 NF 1 BR 1 NR 1 CJE 200 000000E-12 CJC 200 000000E-12 121 CN 2 42 D 87 **** INIȚIAL TRANSIENT SOLUTION TEMPERATURE = 27 000 DEG C NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE ( 1) 5 0000 ( 2) 7521 ( 3) 0 0000 ( 4) 4 7521 ( 5) 9 0000 VOLTAGE SOURCE CURRENTS NAME CURRENT V1 -2 124E-04 V3 0 000E+00 V2 -4 248E-02 TOTAL POWER DISSIPATION 3 83E-01 WATTS **** FOURIER ANALYSIS TEMPERATURE = 27 000 DEG C FOURIER COMPONENTS OF TRANSIENT RESPONSE V(5,4) DC COMPONENT = 2 677038E-01 HARMONIC FREQUENCY FOURIER NORMALIZED PHASE NORMALIZED NO (HZ) COMPONENT COMPONENT (DEG) PHASE (DEG) 1 1 000E+03 5 345E-01 1 000E+00 1 600E-01 0 000E+00 2 2 000E+03 5 316E-01 9 947E-01 -8 959E+01 -8 975E+01 3 3 000E+03 5 269E-01 9 858E-01 -1 794E+02 -1 795E+02 4 4 000E+03 5 203E-01 9 736E-01 9 083E+01 9 067E+01 5 5 000E+03 5 120E-01 9 580E-01 1 041E+00 8 813E-01 6 6 000E+03 5 019E-01 9 392E-01 -8 875E+01 -8 891E+01 7 7 000E+03 4 903E-01 9 173E-01 -1 785E+02 -1 787E+02 8 8 000E+03 4 771E-01 8 927E-01 9 168E+01 9 152E+01 9 9 000E+03 4 625E-01 8 653E-01 1 897E+00 1 737E+00 10 1 000E+04 4 465E-01 8 355E-01 -8 788E+01 -8 804E+01 TOTAL HARMONIC DISTORTION = 2 791763E+02 PERCENT 122 1 0V 0Hz 10KHz 20KHz 30KHz 40KHz 50KHz o V(5,4) Frequency o V(5,4) Frequency Problema 2 **** CIRCUIT DESCRIPTION *etR61 cir CIR cu SWITCH options reltol=5 e-7 itl5=0 numdgt=8 V1 1 0 + DC 100 * V4 4 0 PULSE (0 1 50 u 1 p 1 p 05m 1m ) * L1 1 2 20m C1 3 0 1 25u R1 3 0 400 * d1 2 3 dioda model dioda d 123 * S1 2 0 4 0 S1 * MODEL S1 VSWITCH ( RON=1, ROFF=1MEG, + VON=1, VOFF=0 ) * TRAN 10u 5m * TRAN 10u 5m 4 5m * FOUR 10 K V(3),V(4) PROBE END **** Diode MODEL PARAMETERS dioda IS 10 000000E-15 **** V oltage Controlled Switch MODEL PARAMETERS S1 RON 1 ROFF 1 000000E+06 VON 1 VOFF 0 **** INITIAL TRANSIENT SOLUTION TEMPERATURE = 27 000 DEG C NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE ( 1) 100 00000000 ( 2) 100 00000000 ( 3) 99 20229970 ( 4) 0 00000000 VOLTAGE SOURCE CURRENTS NAME CURRENT V1 -2 481E-01 V4 0 000E+00 TOTAL POWER DISSIPATION 2 48E+01 WATTS 124 110 0 mV 8 0 0mV 4 0 0mV 0V 4 5ms 4 6ms 4 8ms 5 0ms o V(4) Ti me x V(3) Ti me Freque ncy 125 2 0 0V 100V OV-I 1 0Hz o V(3) 1 0 KH z 2 0 KH z 3 0KHz 4 0 KH z Frequency o V(3) F r eque nc y 126 LUCRAREA 7 Problema 1 **** CIRCUIT DESCRIPTION *circ in reg tranz cetrl cir options reltol=1e-6 itl5=0 numdgt=8 NOMOD NOPAGE NOECHO R1 0 1 3 L1 1 2 1 IC=- 10E+01 C1 2 0 1 IC= 10E+01 TRAN 10M 15 SKIPBP PRINT TRAN V(2),I(L1) PROBE END **** TRANSIENT ANALYSIS TEMPERATURE = 27 000 DEG C TIME V(2) I(L1) 1 0000000E-02 2 0000000E-02 3 0000000E-02 4 0000000E-02 5 0000000E-02 6 0000000E-02 7 0000000E-02 8 0000000E-02 9 0000000E-02 1 0000000E-01 1 4960000E+01 1 4970000E+01 1 4980000E+01 1 4990000E+01 1 5000000E+01 9 9011134E-01 9 8041165E-01 9 7091674E-01 9 6166173E-01 9 5251726E-01 9 4349691E-01 9 3484698E-01 9 2619704E-01 9 1754711E-01 9 0942287E-01 2 3808301E-03 2 3716028E-03 2 3623756E-03 2 3531484E-03 2 3439212E-03 9 8027822E-01 9 6102829E-01 9 4228926E-01 9 2414761E-01 9 0628115E-01 8 8872298E-01 8 7208478E-01 8 5544659E-01 8 3880839E-01 8 2347094E-01 9 0939618E-04 9 0587170E-04 -9 0234722E-04 -8 9882274E-04 -8 9529826E-04 127 + I(I1) O V(2) Time + I(L1) V(2 ) **** CIRCUIT DESCRIPTION *circ in reg tranz cetrl cir options reltol=1e-6 itl5=0 numdgt=8 NOMOD NOPAGE NOECHO R1 0 1 3 L1 1 2 1 IC= 10E+01 C1 2 0 1 IC= 10E+01 TRAN 10M 15 SKIPBP PROBE END 128 Time V !(Ы) V(2) **** CIRCUIT DESCRIPTION *circ in reg tranz cetrl cir options reltol=1e-6 itl5=0 numdgt=8 NOMOD NOPAGE NOECHO R1 0 1 3 L1 1 2 1 IC= 10E+01 C1 2 0 1 IC=- 10E+01 129 TRAN 10M 15 SKIPBP PROBE END V I(L1) 0 V(2) Time 0 I(L1) V(2) **** CIRCUIT DESCRIPTION *circ in reg tranz cetrl cir options reltol=1e-6 itl5=0 numdgt=8 NOMOD NOPAGE NOECHO R1 0 1 -3 L1 1 2 1 IC= 10E+01 C1 2 0 1 IC= 10E+01 130 TRAN 1M 2 SKIPBP PROBE END Time V(2 ) 131 V(2 ) **** CIRCUIT DESCRIPTION *circ in reg tranz cetrl cir options reltol=1e-6 itl5=0 numdgt=8 NOMOD NOPAGE NOECHO R1 0 1 3 L1 1 2 1 IC= 10E+01 C1 2 0 -1 IC= 10E+01 TRAN 1M 15 SKIPBP PROBE END Time 132 V(2 ) **** CIRCUIT DESCRIPTION *circ in reg tranz cetrl cir options reltol=1e-6 itl5=0 numdgt=8 NOMOD NOPAGE NOECHO Rl 0 l l L1 1 2 1 IC= 10E+01 C1 2 0 1 IC= 10E+01 TRAN 1M 15 SKIPBP PROBE END Time 133 V(2) **** CIRCUIT DESCRIPTION *circ in reg tranz cetrl cir options reltol=1e-6 itl5=0 numdgt=8 NOMOD NOPAGE NOECHO Rl 0 l -l L1 1 2 1 IC=- 10E+01 C1 2 0 1 IC=- 10E+01 TRAN 1M 15 SKIPBP PROBE Time 134 V(2 ) V(2 ) **** CIRCUIT DESCRIPTION *circ in reg tranz cetrl cir options reltol=1e-6 itl5=0 numdgt=8 NOMOD NOPAGE NOECHO Rl 0 l ln L1 1 2 1 IC= 10E+01 C1 2 0 1 IC= 10E+01 TRAN 1M 15 SKIPBP PROBE 135 Time O I(L1) V( 2 ) Problema 2 **** CIRCUIT DESCRIPTION *circ nelin de ord 2 cnd2 cir L1 0 1 1 IC=6 G1 1 0 TABLE {V(1)}= + (-10,-8),(-1,1),(1,-1),(10,8) TRAN 1M 7 SKIPBP PROBE END 136 2 0A- 0A -2 0A- -4 0V V I(I1) -2 QV QV 2 0V 4 0V 6 0V V(1) *circ nelin de ord 2 cnd2 cir L1 0 1 1 IC=6 G1 1 0 TABLE {V(1)}= + (-10,-8),(-1,1),(1,-1),(10,8) C1 1 0 1M IC=8 TRAN 1M 7 SKIPBP PROBE END V1(C1) 137 *circ nelin de ord 2 cnd2 cir L1 0 1 1 IC=6 G1 1 0 TABLE {V(1)}= + (-10,-8),(-1,1),(1,-1),(10,8) C1 1 0 1M IC=8 TRAN 1M 7 SKIPBP PROBE OPTIONS RELTOL=1e-7 END VL(CL) 138 LUCRAREA 8 Problema 1a **** CIRCUIT DESCRIPTION *el81 cir circ el nelin reg tran options reltol=1e-6 itl5=0 numdgt=8 r1 1 2 15 l1 2 3 10m V1 1 0 sin(0 0 5 1 e+5 0 0 0) d 3 0 dioda model dioda d(Is=8 3e-15 Rs=9 6 + Cjo=300p Vj= 75 + M= 4 Tt=4 e-06) tran 5n 400u 360u * tran 2u 200u 100u 1u probe end **** Diode MODEL PARAMETERS dioda IS 8 300000E-15 RS 9 6 TT 4 000000E-06 CJO 300 000000E-12 VJ 75 M 4 **** INITIAL TRANSIENT SOLUTION TEMPERATURE = NODE VOLTAGE NODE VOLTAGE NODE VOLTAGE ( 1) 0 00000000 ( 2) 0 00000000 ( 3) 0 00000000 VOLTAGE SOURCE CURRENTS NAME CURRENT V1 0 000E+00 27 000 DEG C NODE VOLTAGE 139 400uA 370us 3 9 0us 4 00us 3 6 0us □ I(11) 3 80us Ti me □ I(11) V( 1 ) 140 Problema 1b v=1 5V □ i(li) Ti me □ I(li) V( 1 ) 141 Problema 1c v=6 5V □ i(li) Ti me □ I(li) V( 1 ) 142 Bibliografie M Nifescu, F Constantinescu, ” Bazele electrotehnicii Partea I Teoria circuitelor electrice ” Editura Printech (Andor Tipo), București, 1998, ISBN 973-98453-5-5 F Constantinescu, M Nifescu ” Electrothenică pentru Facultatea de Automatică și Calculatoare”, curs publicat pe internet, M Nifescu, F Constantinescu” Electrothenică pentru Facultatea Energetică”, curs publicat pe internet, A Vladimirescu, ” SPICE ” Editura tehnică București, 1999, ISBN 973-31-1225-9 L O Chua, C A Desoer, E S Kuh, ”Linear and nonlinear circuits” Mc Graw-Hill, 1987 143